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ABSTRACT 

Cold  molecules  and  molecular  ions  are  leading  to  a  renaissance  in  the  field  of  molecular  spectroscopy  just  as  laser- 
cooled  atoms  resulted  in  a  renaissance  of  atomic  physics.  Cold  molecules  enable  the  performance  of  spectroscopy 
with  unprecedented  precision.  Spectroscopy  with  cold  molecules  is  driven  by  many  of  its  modem  ap-  plications  such 
as  precision  measurements,  cold  chemistry  and  quantum  infonnation.  Molecular  ions  trapped  in  RF  Paul  traps  and 
sympathetically-cooled  with  laser-cooled  atomic  ions  have  been  shown  to  be  a  great  platform  to  measure 
spectroscopic  transi-  tions  lines  at  a  precision  beyond  traditional  methods,  where  high  ion  density  is  neces-  sary  but 
difficult  to  achieve  with  the  classical  preparations.  In  this  thesis,  we  perform  vibrational  spectroscopy  on  the  ??  = 
10??=0  and  ??  =  9  ?  ?  =  0  overtone  of  a  trapped  and  sympathetically-cooled  40CaH+  molecular  ion  using  a 
resonance-  enhanced  two-photon  dissociation  scheme.  Our  experiments  are  motivated  by  theo-  retical  work  that 
proposes  comparing  the  vibrational  overtones  of  40CaH+  with  elec-  tronic  transitions  in  atoms  to  detect  possible  time 
variation  of  in  the  mass  ratio  of  the  proton-to-electron.  Due  to  a  lack  of  experimental  data,  we  start  the  search  using  a 
broadband  femtosecond  Ti:Sapph  laser  and  the  guidance  of  theoretical  calculations.  Our  initial  spectroscopy  opens 
the  path  for  more  precise  measurements  to  detect  the  variation  of  proton-to-electron  mass  ratio  and  to  search  for 
astrophysical  presence  of  40CaH+.  Our  method  fits  a  wide  range  of  molecular  ions  that  can  be  co-trapped  and 
efficiently  cooled  by  atomic  ions  and  is  the  first  step  to  higher  precision  mea-  surements.  This  enabled  us  to  measure 
high-order  vibrational  overtones  of  a  molecule  using  only  less  than  300  molecules  for  the  entire  spectrum.  High 
precision  measurements  on  40CaH+  are  enhanced  by  cooling  down  the  motional  temperature  ion  down  to  the 
motional  ground  state.  Future  high  precision  measurements  will  be  conducted  on  narrow  transitions  using  techniques 
that  transforms  the  internal  state  changes  into  changes  in  temperature  of  the  ions  of  40CaH+. 
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SUMMARY 


Cold  molecules  and  molecular  ions  are  leading  to  a  renaissance  in  the  field  of 
molecular  spectroscopy  just  as  laser-cooled  atoms  resulted  in  a  renaissance  of  atomic 
physics.  Cold  molecules  enable  the  performance  of  spectroscopy  with  unprecedented 
precision.  Spectroscopy  with  cold  molecules  is  driven  by  many  of  its  modern  ap¬ 
plications  such  as  precision  measurements,  cold  chemistry  and  quantum  information. 
Molecular  ions  trapped  in  RF  Paul  traps  and  sympathetically-cooled  with  laser-cooled 
atomic  ions  have  been  shown  to  be  a  great  platform  to  measure  spectroscopic  transi¬ 
tions  lines  at  a  precision  beyond  traditional  methods,  where  high  ion  density  is  neces¬ 
sary  but  difficult  to  achieve  with  the  classical  preparations.  In  this  thesis,  we  perform 
vibrational  spectroscopy  on  the  //  =  10^— n  =  0  and  u'  —  9  4—  v  —  0  overtone 
of  a  trapped  and  sympathetically-cooled  40CaH+  molecular  ion  using  a  resonance- 
enhanced  two-photon  dissociation  scheme.  Our  experiments  are  motivated  by  theo¬ 
retical  work  that  proposes  comparing  the  vibrational  overtones  of  40CaH+  with  elec¬ 
tronic  transitions  in  atoms  to  detect  possible  time  variation  of  in  the  mass  ratio  of 
the  proton-to-electron.  Due  to  a  lack  of  experimental  data,  we  start  the  search  using 
a  broadband  femtosecond  Ti:Sapph  laser  and  the  guidance  of  theoretical  calculations. 
Our  initial  spectroscopy  opens  the  path  for  more  precise  measurements  to  detect  the 
variation  of  proton-to-electron  mass  ratio  and  to  search  for  astrophysical  presence 
of  40CaH+.  Our  method  fits  a  wide  range  of  molecular  ions  that  can  be  co-trapped 
and  efficiently  cooled  by  atomic  ions  and  is  the  first  step  to  higher  precision  mea¬ 
surements.  This  enabled  us  to  measure  high-order  vibrational  overtones  of  a  molecule 
using  only  less  than  300  molecules  for  the  entire  spectrum.  High  precision  measure¬ 
ments  on  40CaH+  are  enhanced  by  cooling  down  the  motional  temperature  ion  down 
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to  the  motional  ground  state.  Future  high  precision  measurements  will  be  conducted 
on  narrow  transitions  using  techniques  that  transforms  the  internal  state  changes  into 
changes  in  temperature  of  the  ions  of  40CaH+. 
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CHAPTER  I 


INTRODUCTION 


In  the  field  of  molecular  spectroscopy,  diatomic  spectroscopy  has  yielded  great  divi¬ 
dends  in  our  understanding  of  chemistry  and  physics.  Because  of  the  simple  structure 
of  diatomics,  experimental  manipulation  and  control  is  relatively  easy  compared  to 
complex  molecules.  The  study  of  diatomic  spectra  was  encapsulated  by  Herzberg’s 
seminal  book,  Spectroscopy  of  Diatomic  Molecules  in  1950,  which  summarized  meth¬ 
ods  of  electronic  spectroscopy  measurements  of  diatomics  at  the  time  [1],  The  impact 
of  the  book  has  been  far  reaching,  and  it  is  considered  a  classic. 

Diatomic  molecules  are  attractive  because  relatively  accurate  solutions  of  the 
Schrodinger  equation  can  be  obtained  without  requiring  expensive  computational 
time.  Diatomic  spectroscopy  is  used  as  the  basis  for  the  determination  of  molecular 
structure  and  reaction  dynamics.  The  simple  nature  of  diatomic  molecules  allows  for 
high  precision  ab  inititio  calculations  because  there  are  few  electrons  to  worry  about 
versus  complex  molecules.  Improvements  on  experimental  determination  of  molecular 
constants  serves  as  a  benchmark  on  which  molecular  theory  can  be  refined. 

Not  only  is  the  study  of  molecular  spectra  important  in  determining  molecular 
structure,  but  it  became  apparent  that  diatomic  molecules  play  an  important  role  in 
the  chemical  and  physical  processes  occurring  in  the  interstellar  medium.  Diatomic 
molecules  serve  as  the  building  blocks  for  larger  molecules  in  the  ISM.  An  example  is 
the  formation  of  the  Hjj"  molecule,  which  is  believed  to  be  formed  by  reactive  collisions 
between  and  H  [2],  Recently,  it  has  been  suggested  that  some  of  the  fundamental 
constants  of  nature  such  as  the  fine  structure  constant,  a,  and  the  proton-to-electron 
mass  ratio,  p,  may  not  be  constant  over  cosmological  time  scales  [3].  High  precision 


1 


measurements  on  molecular  spectra,  can  test  the  temporal  variation  of  fundamental 
constants  of  nature.  Diatomic  molecules  are  posed  as  ideal  candidates  for  the  study 
of  fundamental  physics. 

1.1  Diatomics  in  interstellar  and  precision  measurements 

Extraterrestrial  spectroscopy  traces  its  origins  to  German  optician  and  physicist, 
Joseph  von  Fraunhofer,  who  constructed  high  quality  lenses  to  make  telescopes,  which 
he  used  to  observe  dark  lines  of  the  solar  spectrum  in  1718  [4].  He  invented  a  he¬ 
liometer  and  a  diffraction  grating  and  studied  lines  previously  noticed  by  W.H.  Wol¬ 
laston  [5].  Years  later,  in  1769,  G.  R.  Kirchhoff  and  R.  W.  Bunsen  demonstrated 
that  the  comparison  of  laboratory  spectra  with  emission  spectra  from  stars  was  a 
useful  tool  for  studying  the  cosmos  [6].  This  signified  the  power  of  spectroscopy  as 
an  instrument  to  study  chemical  and  physical  properties  of  interstellar  processes  and 
fundamental  physics.  The  technique  of  matching  laboratory  spectra  with  absorption 
spectra  from  interstellar  medium  (ISM)  is  still  being  used  today  to  characterize  many 
lines  from  the  ISM  [7,  8].  While  significant  progress  has  been  made  to  identify  many 
emission  lines,  quite  a  number  still  remain  uncharacterized,  often  referred  to  as  “in¬ 
terstellar  weeds” .  This  presents  great  opportunities  to  study  molecular  spectra  with 
applications  to  astrochemistry  and  astrophysics. 

Research  in  interstellar  chemistry  started  when  the  first  molecules  were  detected 
in  space  in  the  1930s  [9].  The  first  one  was  the  CH  radical  and  sometime  later  the 
CH+  ion  and  CN  molecule  [10,  11].  These  had  an  electronic  spectra  in  the  region  of 
the  atmosphere  that  was  transparent  to  the  earth.  The  spectra  from  diffuse  clouds 
is  dominated  by  atomic  species  and  diatomic  molecules.  Later  discoveries  of  complex 
molecules  indicated  that  diatomic  molecules  played  an  important  role  in  the  formation 
of  those  molecules.  Since  atomic  hydrogen  is  the  most  abundant  atom  in  space  (~ 
70%),  it  is  no  surprise  that  many  molecules  in  the  ISM  are  hydrides. 
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Hydride  molecular  ions  are  of  great  interest  because  many  of  their  neutral  coun¬ 
terparts  have  been  discovered  in  regions  that  favor  the  formation  of  ions.  Several 
hydride  molecular  ions  have  been  detected  in  the  solar  photospheric  spectrum  includ¬ 
ing  MgH+,  OH+,  SiH+  and  A1H+  [12,  13].  These  molecules  are  also  present  in  comets, 
cool  stellar  atmospheres,  and  interstellar  media.  Due  to  the  large  ionization  potential 
of  these  atoms,  it  is  believed  that  these  hydrides  will  be  formed  deep  inside  the  solar 
atmosphere.  These  molecules  become  the  basis  of  weak  unresolved  lines  with  enhanced 
opacities.  Studies  of  collisional  processes  involving  highly  excited  states  of  atomic  ions 
such  as  Ca+  and  Mg+  with  H2  in  stellar  atmospheres  leads  to  the  broadening  of  the 
spectra  of  these  atomic  ions  in  the  resonance  doublet  solar  spectrum.  [14,  15].  The 
need  for  more  investigations  of  molecular  ions  in  the  ISM  was  called  for  by  Herzberg 
in  his  essay  in  the  book  Molecular  Ions:  Spectroscopy,  Structure,  and  Chemistry  when 
he  said: 

Although  the  number  of  molecular  ions  observed  in  the  interstellar  medium 
is  still  small  it  is  generally  agreed,  on  the  basis  of  the  work  by  Herbst  and 
Klemperer,  Watson  and  others,  that  ions  play  a  major  role  in  the  chemi¬ 
cal  reactions  that  take  place  in  the  interstellar  medium.  This  is  so  largely 
because  ion-molecule  reactions  are  fast  reactions;  they  have  no  activation 
energy. 

The  above  statement  highlights  the  general  importance  of  molecular  ions  in  driving 
the  processes  in  the  ISM  while  also  lamenting  on  the  lack  of  observed  spectra  [16]. 

Neutral  CaH  has  been  observed  in  sunspots  where  the  temperatures  are  above 
the  ionization  potential  of  /( Ca)=  6.113  eV  [17].  Furthermore,  the  observation  of 
MgH+  lines  in  the  solar  spectrum  while  Mg  has  an  ionization  potential  of  /( Mg)  = 
7.646  eV  indicates  that  40CaH+  should  be  observed  in  the  solar  spectrum  [18].  Ca+ 
has  been  observed  in  the  atmosphere  where  it  is  believed  to  be  deposited  by  meteor 
ablation  [19].  However,  despite  several  calls  for  the  spectra  for  40CaH+,  there  had  not 
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been  any  experimental  spectra  of  40CaH+.  Few  theoretical  studies  of  40CaH+  exist, 
which  form  the  basis  of  this  experiment  [17,  20-25]. 

1.2  Diatomics  and  fundamental  physics 

Temporal  variations  of  fundamental  constants  such  as  a,  the  fine  structure  constant, 
and  /i,  the  proton-to-electron  mass  ratio,  are  important  in  testing  current  physical 
theories  and  development  of  new  theories  such  as  supersymmetry.  The  local  position 
invariance  principle  of  the  Einstein  equivalence  principle  states  that  there  should  be 
no  variance  of  fundamental  constants  of  nature  in  time  or  position  [26,  27].  How¬ 
ever,  proposed  theories  allow  for  the  variation  in  space  and  time  of  fundamental 
constants  [28] .  These  models  deriving  from  string  theory  include  space-time  with  ex¬ 
tra  dimensions  where  the  geometry  can  vary  as  well  as  the  the  amplitude  of  light 
scalar  fields  that  can  couple  to  ordinary  matter.  These  variations  change  the  ap¬ 
parent  values  of  constants  [29].  These  models  are  important  in  the  bid  to  explain 
the  dark  energy  in  the  universe.  Testing  these  theories  requires  precision  measure¬ 
ments  of  atomic  and  molecular  spectra.  Diatomic  molecules  have  been  proposed  as 
ideal  candidates  to  study  the  possible  time  variation  of  the  proton-to-electron  mass 
ratio,  y.  Prior  to  high  precision  spectroscopy,  the  time  variation  of  fundamental  con¬ 
stants  was  determined  by  comparing  red-shifted  spectra  from  quasars  and  current 
laboratory  spectra  [30].  However,  such  measurements  only  have  an  accuracy  of  1CT6, 
which  is  significantly  lower  than  the  needed  accuracy  for  the  detection  of  year-to-year 
variation. 

Significant  progress  has  been  made  in  atomic  spectroscopy  to  place  an  upper  bound 
on  the  variation  of  a.  Rosenband  et  al.  compared  the  ratio  of  Al+  to  Hg+  electronic 
transitions  to  an  accuracy  of  5.2xl0~17,  and  the  change  in  a  to  be  (1.6±2.3)  x  10-1' 
per  year  [31].  The  time  variation  in  y  cannot  be  easily  detected  in  atomic  ions  because 
y  weakly  depends  on  isotopic  effects  in  atomic  ions.  On  the  other  hand,  the  time 
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variation  in  /i  scales  with  the  transitions  in  molecules  as  Eei  :  Evib  :  Erot  «  1  :  yU~1//2  : 
/2~1.  Comparing  the  rovibrational  spectra  of  molecules  can  offer  great  insight  on  the 
time  variation  of  //.  There  have  been  many  proposals  and  experiments  to  determine 
the  upper  bound  of  //  [28,  29,  32,  33].  So  far  none  of  these  experiments  have  been 
able  to  achieve  the  required  precision  of  greater  than  10-15  needed  to  detect  the 
variation  on  a  year-to-year  basis.  A  recent  experiment  has  been  able  to  measure  the 
mass  of  the  electron  by  comparing  the  cyclotron  frequency  of  a  heavy-ion  system 
to  the  frequency  of  the  precision  of  the  electron  spin  [34],  However,  comparing  the 
rovibrational  spectra  of  molecules  remains  the  best  method  to  detect  yearly  variations 
in  //.  For  example,  40CaH+  has  been  proposed  to  measure  the  temporal  variation  in 
fundamental  constants  [33]. 

1.3  Cold  molecular  ions 

Despite  the  fact  that  many  other  molecular  ions  have  been  identified  in  the  interstellar 
medium,  research  in  molecular  spectroscopy  has  largely  favored  neutral  molecules 
over  molecular  ions.  This  is  because  traditional  spectroscopy  techniques  require  large 
molecular  densities,  ^108  molecules/cm3  to  obtain  a  good  signal-to-noise  ratio.  This 
becomes  difficult  with  ions  because  Coulombic  repulsion  between  the  ions  limits  ion 
density  [35].  Furthermore,  detection  schemes  used  in  neutral  molecule  studies  requires 
the  production  of  an  ion,  which  can  be  detected  as  a  result  of  probing  a  transition,  are 
not  very  useful  for  molecular  ions.  Consequently,  there  are  many  unassigned  lines  in 
the  interstellar  spectrum  [36].  Diatomic  molecular  ions  could  provide  some  solutions 
to  the  observed  but  unaccounted  interstellar  lines.  To  get  around  these  limitations,  ion 
trapping  and  Doppler  cooling  can  perform  high  precision  measurements  on  molecular 
ions,  because  it  removes  line  broadening  processes  such  as  pressure  and  Doppler 
broadening. 

Cold  molecular  ions  have  an  advantage  over  neutral  molecules  because  they  can  be 
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stored  with  a  combination  of  RF  and  DC  fields.  This  allows  for  longer  interrogation 
times,  Doppler  free,  and  collision  free  measurements.  In  addition,  by  trapping  ions  as 
a  linear  crystal  chain  at  the  center  of  the  RF  null,  the  Stark  shift  can  be  reduced  so 
that  measurements  can  be  done  with  an  accuracy  greater  than  10-15  [29]. 

The  combination  of  radio  frequency  traps  and  laser  cooling  seeks  to  overcome 
the  challenges  of  measuring  the  spectra  of  molecular  ions.  In  ion  traps,  the  ions 
are  confined  in  a  specified  region  for  a  long  time  where  lasers  can  address  them 
repeatedly.  Thus,  the  lack  of  density  is  no  longer  a  problem  since  the  position  of 
the  ion  is  exactly  known,  therefore,  the  density  becomes  effectively  infinite  [37].  In 
addition,  removing  the  translational  energy  of  the  ion,  i.e.  laser  cooling,  removes  the 
excess  energy,  which  allows  for  weak  transitions  to  be  probed  with  a  high  signal-to- 
noise  ratio.  Because  the  ions  are  trapped  in  isolated  environments  with  pressures  lower 
than  10~9  torr,  the  collision  and  pressure  broadening  effects  are  avoided.  The  40CaH+ 
molecular  ion  could  help  characterize  some  of  these  interstellar  weeds.  While  most 
diatomic  molecular  ions  cannot  be  directly  laser-cooled,  they  can  be  cooled  through 
their  Coulombic  interaction  with  laser-cooled  atomic  ions  in  a  trap,  a  process  called 
sympathetic  cooling  [38]. 

1.4  Sympathetically-cooled  molecular  ions 

The  advent  of  laser  cooling  produced  many  possibilities  for  scientific  research.  How¬ 
ever,  not  many  ions  can  be  directly  laser-cooled  because  laser  cooling  requires  closed 
optical  cycles  that  can  be  pumped  repeatedly  without  decay  to  other  metastable 
states1.  Other  ions  have,  very  deep  UV  regions  for  optical  cycling,  which  becomes  ex¬ 
perimentally  unfeasible.  Many  molecules  have  a  similar  fate;  due  to  vibrational  and  ro¬ 
tational  structure  of  molecules,  direct  laser  cooling  is  not  possible  for  most  molecules. 

1A  recent  paper  [39]  suggests  that  cooling  can  be  achieved  without  spontaneous  emission  by 
restricting  the  atom-light  interaction  to  a  time  short  compared  to  the  cycle  of  absorption  followed 
by  natural  decay.  The  applicability  of  this  method  in  trapped  ions  has  not  yet  been  demonstrated. 
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In  order  to  directly  cool  molecules,  the  internal  structure  of  the  molecule  needs  strong 
Franck-Condon  overlap  to  minimize  decays  to  other  rovibrational  levels.  This  means 
that  a  molecule  can  be  laser-cooled  with  few  lasers.  Neutral  molecules  such  as  SrF, 
YO,  and  CaF,  have  diagonal  Franck-Condon  overlap,  which  allows  them  to  be  laser- 
cooled  [40-43].  For  molecular  ions,  only  a  few  candidates  have  similar  Franck-Condon 
overlap  to  SrF  that  have  been  identified,  A1H+,  BH+  and  SrO+  [44,  45]. 

For  most  molecules,  sympathetic  cooling  with  a  laser-cooled  atomic  ion  remains 
the  solution  for  reducing  the  motional  energy  and  performing  precision  spectroscopy 
[46].  Sympathetic  cooling  involves  multiple  Coulomb  interacting  ions  stored  in  trap 
where  one  ion  species  is  laser-cooled  and  the  other  is  not  laser-cooled.  The  laser- 
cooled  ions  bring  all  the  ions  in  the  trap  to  thermal  equilibrium  with  the  temperature 
of  the  laser-cooled  ion  through  their  mutual  Coulombic  interaction  [47-49].  The 
laser-cooled  and  the  sympathetically-cooled  ion  pair  form  ordered  structures  known 
as  Coulomb  crystals  [50].  Sympathetic  cooling  does  not  disturb  the  internal  structure 
of  the  molecular  ion  because  the  long-range  Coulomb  interaction  energy  between 
the  trapped  ions  is  far  less  than  rovibrational  energy  spacing.  Thus,  the  internal 
energy  of  trapped  molecular  ions  is  equilibrated  to  background  black  body  radiation 
(BBR),  which  is  at  room  temperature.  Because  there  is  no  need  for  direct  laser¬ 
cooling,  sympathetic  cooling  increases  the  catalog  of  ions  that  can  be  cooled  and 
trapped.  Generally,  sympathetically-cooled  ions  and  laser-cooled  ions  do  not  have 
similar  transitions  such  that  precision  measurements  on  the  sympathetically-cooled 
ions  do  not  disturb  the  laser  cooling  of  the  atomic  ions  [51,  52], 

Many  studies  have  been  conducted  with  sympathetically-cooled  molecular  ions 
including  isotope  effects  in  the  reaction  between  Mg+  and  HD  [53],  charge  exchange 
reactions  [54],  velocity  selected  collisions  [55,  56],  sympathetically-cooled  biomolecules 
[57-60],  rovibrational  spectroscopy  [61-65],  and  internal  cooling  of  molecules  [66 — 
73] .  Sympathetically-cooled  ions  are  also  applicable  to  quantum  information  because 
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sympathetically-cooled  ions  are  rarely  perturbed  by  the  Doppler  cooling  process  of 
the  laser-cooled  ions  [74-77].  For  the  vibrational  spectroscopy  of  40CaH+,  we  use 
sympathetically-cooled  40CaH+  ions  with  laser-cooled  40Ca+  in  a  Paul  trap. 

Despite  the  difficulty  in  experimenting  with  molecular  ions,  several  techniques 
have  been  employed  over  the  years  to  measure  molecular  ion  spectra.  Some  tradi¬ 
tional  spectroscopy  techniques  such  as  crossed  beams  were  implemented  but  with 
limited  success  [16].  Ion  fragmentation  or  photodissociation  has  received  wide  spread 
applicability  in  molecular  ion  spectroscopy  because  of  the  ease  to  which  the  technique 
can  be  performed  [78].  In  ion  fragmentation  spectroscopy,  one  or  more  lasers,  tuned 
to  a  state  of  interest,  irradiates  the  ion.  If  the  photons  are  in  resonance  with  the 
transition  of  interest,  the  molecule  will  absorb  the  photons  and  be  excited  beyond 
its  dissociation  threshold.  The  fragments  can  be  collected  and  analyzed  for  angu¬ 
lar  distribution,  velocity,  and  types  of  fragments,  to  probe  dynamics.  Alternatively, 
the  number  of  fragments  or  speed  of  fragmentation  can  be  probed  as  a  function  of 
the  laser  frequency,  probing  the  rovibrational  spectra  [79].  Several  pioneering  studies 
have  been  conducted  using  this  technique.  These  include  infrared  multiphoton  disso¬ 
ciation  (IR.MPD)  of  polyatomic  ions  [80,  81],  resonant  ion-clip  infrared  spectroscopy 
(RIDIRS)  of  water  complexes  [82,  83],  vibrationally  and  rotationally  resolved  elec¬ 
tronic  spectra  of  metal-ion  complexes,  vanadium  oxide-carbonyl  cations  [19,  84],  and 
IR  clip-ultraviolet  photodissociation  of  protonated  tyrosine  and  phenylalanine  [85].  In 
this  thesis,  we  use  ion-fragmentation  spectroscopy  of  single  trapped  molecular  ions  to 
study  the  vibrational  overtones  of  40CaH+. 

1 . 5  Thesis  organization 

This  thesis  will  be  focused  on  the  vibrational  spectroscopy  of  40CaH+  sympathetically 
laser-cooled  with  40Ca+  ions.  This  is  the  first  such  measurement  of  bound  transitions 
of  40CaH+.  The  rest  of  the  thesis  is  organized  as  follows:  Chapter  2  discusses  the 


theory  of  trapping  and  laser  cooling  ion  in  RF  traps.  Chapter  3  discusses  the  molecu¬ 
lar  orbital  theory  of  40CaH+  found  in  the  literature  including  the  reaction  mechanism 
by  which  the  molecular  ion  is  formed.  Chapter  4  shows  the  kind  of  instruments 
used  to  study  the  ion.  Chapter  5  demonstrates  the  use  of  resonance  enhanced  mul¬ 
tiphoton  dissociation  technique  to  record  for  the  first  time  vibrational  overtones  of 
40CaH+.  Chapter  6  shows  a  proposal  of  an  experiment  where  the  spectra  of  40CaH+ 
is  measured  by  looking  at  the  heating  of  the  sidebands  of  a  ground  stated  cooled 
molecular  ion.  In  the  same  chapter  we  chronicle  the  first  step  towards  precision  mea¬ 
surements  of  by  cooling  the  40CaH+  ion  motion  inside  the  trap  to  the  ground  state 
and  gives  indications  of  future  research  directions  in  the  Brown  lab. 
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CHAPTER  II 


THEORY  ION  TRAPPING  AND  LASER  COOLING 


Key  to  the  experiments  described  in  this  thesis  is  the  trapping  and  cooling  of  charged 
particles.  This  is  facilitated  by  a  Paul  RF  quadrupole  trap  that  confines  the  ions  in 
three  dimensions  using  a  combination  of  RF  and  DC  fields.  Doppler  cooling  is  used 
to  reduce  the  translational  energy  down  to  millikelvin  temperatures.  In  this  chapter 
the  operation  principle  of  a  Paul  trap  and  the  interactions  of  the  ions  with  the  lasers 
is  discussed.  The  reader  is  referred  to  the  following  references  for  more  details  [86-94] 

2.1  Theory  of  trapping  ions 

Trapped  ion  experiments  use  an  RF  Paul  trap  named  after  Wolfgang  Paul  [95].  To 
achieve  trapping  a  quadrupole  held  is  used  to  confine  the  charged  particles  in  space 

$  =  +  Py2  +  7^2),  (1) 

zr0 

where  r0  is  the  radius  of  the  trap  and  a,  (3, 7  are  trapping  coefficients.  If  the  potential 
is  static,  the  point  r0  =  0  is  unstable  because  the  Laplace  condition  A$=0  imposes 
the  requirement  that  a  +  (3  +  7  =  0.  This  results  in  confinement  in  two  of  the 
directions,  but  anti-confinement  in  the  other  because  7  =  —(/?  +  a).  To  overcome 
this,  the  voltage  is  made  periodic  in  time.  An  oscillating  voltage  results  in  focusing 
and  defocussing  in  the  x  and  z  directions.  The  potential  is  applied  with  a  DC  voltage 
U0  and  RF  voltage  V0  at  a  frequency  D 

$0  =  UQ  +  V0  cos  f it.  (2) 

A  schematic  of  the  a  typical  linear  Paul  trap  is  shown  in  Fig.  1. 
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(a)  3D  sketch 


(b)  Cross  section 


Figure  1:  A  schematic  of  a  linear  Paul  trap  configuration,  (a)  shows  a  3D  cartoon  of 
the  trap  where  an  RF  voltage  V0  cos  Qt  is  applied  on  two  diagonally  opposed  electrodes 
while  the  other  two  are  kept  at  RF  ground.  To  trap  along  the  z  axis,  U0  static 
voltage  is  applied  along  near  the  end  of  the  segments.  The  red  segments  represents 
RF  common  segments  which  can  be  biased  with  DC  to  ensure  that  the  ion  sits  on 
the  RF  null,  (b)  shows  a  cross  section  of  trap  indicating  the  radius  of  the  trap. 


The  applied  voltage  on  the  electrodes  gives  the  potential  in  the  trap  described  by, 

kUq  2  x2  +  y2 


x2  —  v2 

$  =  — — Vo  COS  Qt  + 


2'ro 


[z  - 


-) 


(3) 


where  k  is  the  constant  geometric  factor. 

An  illustration  of  the  operation  of  an  oscillating  held  to  trap  an  ion  is  shown  in 
Fig.  2  which  shows  the  polarity  of  the  held  switching  and  the  ion  inside  the  held.  As 
the  ion  moves  in  one  direction,  the  held  polarity  is  switched,  resulting  in  repulsion  of 
the  ion.  When  the  ion  changes  direction,  the  polarity  is  hipped  again  causing  it  to  be 
repulsed  again.  This  makes  the  ion  oscillate  back  and  forth  around  the  center  of  the 
trap.  The  stability  of  the  ion  inside  the  trap  requires  that  ion  moves  slow  relative  to 
the  switching. 
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(a)  -V  (b)  V 


Figure  2:  An  illustration  of  the  oscillating  field  to  trap  ions  showing  the  two  states  as 
the  fields  oscillates  between  positive  and  negative.  The  polarity  of  the  field  is  switched 
continuously  to  cause  the  ion  to  oscillate  about  the  center  of  the  trap.  The  ion  must 
move  slower  than  the  switching  in  order  for  trapping  to  be  realized. 

2.2  Motion  inside  a  trap 

The  equations  of  motion  of  an  ion  inside  the  trap  can  be  a  described  by  a  charge  in 
a  field  governed  by, 

m~^(x,y,z)  =  (4) 

This  is  done  by  explicitly  expressing  each  component  along  the  coordinates, 


e 

x  H - 


771 


y  + 


e 

m 


-^"7  +  Vq a  cos  Qt  )  x  =  0, 


z - (Uorfz  =  0, 

m 


y7  -  V0acosttt  }y  =  0, 


(5) 

(6) 
(7) 
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where  m  is  the  mass  of  the  trapped  particle.  The  ion  feels  a  small  averaging  force 
towards  the  center  because  the  held  is  inhomogeneous.  These  equations  are  the  well 
known  Mathieu  equations,  which  in  dimensional  quantities  take  the  form, 


d2x 

dr2 


(ax  +  2 qx  cos  2 r)x  =  0, 


d2z 

dr2 


(az  —  2 qz  cos  2 r)z  =  0, 


d2y 

dr2 


+  ( ay 


2 qy  cos  2 r)y  =  0, 


where 


4  KeU0  az 


a*-ay~  mrMy2  ~  2  ’ 


qx  =  ~qy 


2eVo 
mrfcY2  ’ 


f It 

Qz  =  0,  T  =  — . 


(8) 

(9) 

(10) 


(11) 

(12) 

(13) 


The  equations  have  stable  regions  and  unstable  regions,  which  are  dependent  on  the 
values  of  a  and  q  but  not  the  velocity  of  the  ion.  In  experiments,  the  ions  are  typically 
confined  in  a  fixed  ratio  of  a/q,  depending  on  the  mass  of  the  ion  and  the  applied 
voltages.  It  is  common  to  show  the  stability  diagram  near  a  =  0  because  most  traps 
are  designed  to  operate  in  this  region.  The  analytical  solution  of  the  Mathieu  equation 
can  be  expressed  as, 

n= oo  n=oo 

u(t )  =  ole?T  C2ne2inT  +  a"e~^  C2ne~2inT,  (14) 

n=— oo  n=— oo 

where  fi  =  a  ±  i/3,  a'  and  a"  are  integration  constants  that  depend  on  the  initial 
conditions,  and  the  coefficients  C2n  satisfy  the  recursion  relation  and  depend  on  the 
values  of  a  and  q.  In  our  experiments,  we  set  the  a—  0,  so  that  for  the  40Ca+  ion,  with 
mass  40  amu,  we  trap  when  0 <q  <0.9  as  shown  Fig.  3. 
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Figure  3:  The  stability  diagram  of  the  quadrupole  field  showing  the  lowest  stability 
region.  Along  the  a=0  axis,  the  ion  is  stable  for  0<g<0.9.  The  Mathieu  equations 
have  a  stable  solution  inside  the  darker  shaded  region.  In  order  to  trap  ions,  Uo  and 
Vo  parameters  are  chosen  carefully  so  that  they  fall  inside  the  stable  part  of  the  a-q 
space. 

The  stability  of  the  equation  is  dependent  on  the  value  of  f3.  If  /3  is  not  an  integer 
and  /d  is  purely  imaginary,  then  the  solutions  of  the  equation  are  stable.  However,  if 
/3  is  an  integer,  the  solutions  are  periodic  but  unstable.  To  solve  the  equation,  we  use 
the  trigonometric  relation, 

el6  =  cos  6  +  i  sin#.  (15) 

If  we  substitute  equation  15  in  to  equation  14,  we  obtain 

n= oo  n= oo 

u{t)  =  aq  C2n  cos(2n  ±  0)t  +  a"  C2n  sin(2n  ±  /3)r,  (16) 

n=— oo  n=— oo 

where  C2n  represents  the  amplitude  of  the  oscillatory  motion  and  (2 n  ±  f3)r  is  the 
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frequency.  If  we  write 


u>nt  =  (2 n  ±  9)t, 
T  =  fit/ 2, 


(17) 

(18) 


we  find  the  following, 


(2 n  ±/3)fl 


P  ~  \la  +  2’ 

The  value  n  —  0  describes  the  fundamental  frequency  of  the  ion  motion; 


(19) 


(20) 


L Oq  —  Pfl/ 2, 


(21) 


substituting  the  value  of  (3,  we  find, 


eVn 


u0  = 


(22) 


\[2mr/p? 

when  o=0.  The  oscillation  of  the  ion  inside  the  trap  can  be  described  by  the  low  fre¬ 
quency  harmonic  oscillation  called  secular  motion  and  the  high  frequency  oscillation 
called  micromotion.  The  potential  describing  the  secular  motion  is  called  the  pseudo¬ 
potential  to  distinguish  it  from  the  trapping  potential.  The  depth  of  this  potential, 
D,  is  described  by, 


m 


D  =  — wr n. 


(23) 


For  3D  ion  traps  described  above,  trap  depths  between  5  eV  to  20  eV  are  easily 
achieved. 


2.3  Atom  interaction  with  radiation 

The  interaction  of  an  atomic  system  with  electromagnetic  radiation  near  resonance 
with  an  atomic  transition  can  be  described  by  the  time  dependent  Schrodinger  equa¬ 
tion, 

ih—  =  H^>.  (24) 
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The  Hamiltonian  describing  a  trapped  atom  in  a  harmonic  well  interacting  with  a 
laser  can  be  expressed  as, 


H  =  H0  +  Hu 


(25) 


where  H0  describes  the  particle  without  the  laser  interaction  and  the  H\  is  the  pertur¬ 
bation  of  the  system  caused  by  the  atom  interaction  with  the  laser.  For  laser  cooling, 
two  levels  are  typically  used  and  hence  will  be  considered  for  the  treatment  of  the 
Doppler  cooling  process.  Suppose  the  system  is  described  by  a  simple  two-level  sys¬ 
tem  with  Hq  describing  the  system  with  wavefunctions  To  and  T  i  and  energies  E0 
and  Ei  respectively  as  in  Fig.  4.  The  energy  gap  is  E\  —  E0  =  hv i0  =  fauio  =  —fvjJoi. 


I - % 


Figure  4:  A  two  level  atomic  system  with  the  wavefunctions  and  the  energy  gap 
between  the  levels.  E0  is  the  energy  of  the  ground  state  T0  and  E\  is  the  energy  of 
the  excited  state  T]. 

The  solutions  to  the  unperturbed  Hamiltonian  take  the  form, 


T0  =  il;0e-iEot/h  =  ip 0e-iuJo\ 


(26) 


q>1  =  ipxe~iExtlh  =  ipxe-^x\ 


(27) 
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where  ay  =  Ei/K.  The  time  dependent  wavefunction  describing  the  perturbed  Hamil¬ 
tonian  becomes  a  linear  combination  of  the  unperturbed  wavefunctions  as, 

*(t)  =  Co^e-^  +  (28) 

where  the  time  dependent  coefficients  Co(t)  and  C\  (t)  satisfy  the  normalization  condi¬ 
tion  (T(f)|T(f))  =  |C*i(f)|2  +  |C2(t)|2  =  1.  The  term  |Cj(t)|2  describes  the  probability 
that  the  atom  will  be  in  the  ith  state.  To  determine  the  evolution  of  the  system,  we 
solve  the  Schrodinger  equation  by  substituting  equation  28  into  equation  24  to  find 

zh[C0(t)^0e-iuot  +  =  H&oiWoe-***  +  (29) 

where  C\  =  Ct/St  is  the  time  derivative  of  the  coefficient.  To  solve  the  Schrodinger 
equation,  we  multiply  by  or  ip^elWlt  and  integrate  over  all  space  to  obtain, 

ihC0  =  C'o(f)('0o|i?i|'0o)  +  C^olHMe-™0*,  (30) 

ihc1  =  +  C0(^i|^i|V’o)eiwlot.  (31) 

Since  the  term  H i  describes  the  dipole  transition  /i  in  the  presence  of  an  electric  held 
eocos(o;f),  making  it  — /reo  cos  (cot),  the  function  is  odd,  where  eo  is  the  amplitude  of 
the  electric  held.  However,  the  products  of  \ipi\2  are  even  functions.  Because  atomic 
states  have  parity  with  the  respect  of  the  space-fixed  coordinates,  (?/y|lhi|^/y)  =  0. 
Thus  the  equations  30  and  31  become, 


ihC0  =  —Ci(t)Mme0e  tU}lot  cos  (cot), 


(32) 


ihC\  =  —  Co(t)Moie0e^olt  cos(o;f). 


(33) 


The  function  Mm  is  the  transition  dipole  moment  given  by  ('01|/u|'j/’o)-  If  we  define  a 
parameter  we  call  the  Rabi  frequency  as, 

MWE 


fin  — 


h  # 


(34) 
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and  then  we  use  the  trigonometry  identity  that  cos  (cat)  =  ~te  ,  we  can  rewrite 
the  equations  32  and  33  as, 


Cn  = 


C\  = 


iCofloie^10-^  +  ei(wi0+t^) 


The  terms  with  e*huo-^)i  and  e-i(aji0-a;)t  represent  the  slow  oscillations  of  the  function 
in  time  while  the  others,  and  e-*(wio +w)*,  represent  fast  oscillations  in  time. 

Because  the  fast  oscillations  average  to  zero,  they  can  be  ignored,  a  process  called 
the  rotating  wave  approximation  (RWA).  If  we  define  the  detuning  as  A  =  u  — 
then  the  solution  to  the  coefficients  is, 


iQ0elACi(t ) 
2 


(37) 


iQoetACo(t) 
~  2  ' 


(38) 


These  equations  can  be  solved  analytically  using  the  initial  conditions  and  we  hnd 
that, 


Cl(t)  =  I(u:)si,1(?)e"A‘/2’  (40) 

where  =  [(Oo)2  +  A2]1//2.  The  probability  that  the  atom  will  be  found  in  the  excited 
state  is  given  by, 

iClWr  =  |^(f).  (4i) 

The  state  of  the  atoms  evolves  in  time  depending  on  the  detuning  of  the  laser  from 
the  transition.  The  on- resonant  evolution  is  shown  in  Fig.  5. 
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Figure  5:  Rabi  oscillations.  The  red  solid  curve  shows  the  probability  of  finding  the 
atom  in  the  excited  state  in  a  two  level  system  driven  by  a  laser  on  resonance  with 
the  atomic  transition  if  there  is  no  damping  of  the  system.  The  blue  dashed  curve 
shows  the  decay  of  coherence  between  the  two  states. 

2.4  Doppler  cooling 

The  removal  of  translational  temperature  of  the  ion  requires  optical  cycling  on  an 
ion  transition  with  a  short  lifetime,  typically  <100  ns  for  practical  experimental 
purposes.  This  is  achieved  by  using  broad,  fast  decaying  transitions  that  will  al¬ 
low  fast  cooling  times.  While  laser  cooling  can  be  achieved  on  narrow  metastable 
transitions,  these  transitions  cycle  photons  slowly  to  be  useful  for  the  first  stage  of 
cooling  [96,  97].  Laser  cooling  takes  advantage  of  the  Doppler  effect.  A  red-detuned 
laser  traveling  in  the  opposite  direction  of  an  ion  is  absorbed  by  the  ion  because  the 
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frequency  is  Doppler  shifted  to  be  on  resonance  with  the  atomic  transition  [98].  The 
absorbed  photon  is  spontaneously  emitted  by  the  ion  resulting  in  a  momentum  kick, 
photon  recoil.  After  many  photons  (psIO6)  have  been  scattered,  the  momentum  from 
the  photon  recoil  is  averaged  out  resulting  in  a  net  slowing  down  of  the  ion,  i.e., 
translational  cooling.  This  proceeds  until  the  heating  from  the  random  momentum 
kick  and  the  cooling  process  cancel  each  other  out  resulting  the  Doppler  cooling  limit. 

If  the  laser  is  slightly  detuned  from  the  atomic  transition  such  that  u  <  0Jqi  with  a 
wave  vector  k,  and  it  moves  with  a  velocity  v  in  the  direction  opposite  to  the  direction 
of  travel  of  the  atom,  the  Doppler  shifted  frequency  of  the  photon  becomes  resonant 
with  the  transition  through  (ta  — k.v)  =  caoi-  Since  a  short  lived  excited  state  is  chosen 
for  Doppler  cooling,  the  excited  state  will  spontaneously  emit  photons  in  a  random 
direction.  The  radiation  field  will  gain  an  energy  of  ft|k.v|  while  the  atom  will  lose  the 
same  amount  in  the  form  of  kinetic  energy  due  to  the  energy  conservation.  Repeating 
this  cycle  many  times  results  in  lower  kinetic  energy  and  net  cooling  of  the  ion  due 
to  low  kinetic  energy. 

The  velocity  capture  range  Av  of  the  Doppler  cooling  process  is  limited  by  the 
natural  linewidth,  T,  of  the  excited  state, 


kAr;  pa  T. 


(42) 


Since  the  radiation  emitted  by  the  atom  is  spontaneous  and  results  in  heating  and 
cooling  of  the  atom,  the  achievable  temperature  of  Doppler  cooling  is  limited  by  the 
competition  between  the  cooling  and  the  heating  processes  [99].  The  lowest  achievable 
temperature  becomes, 


Td  = 


hr 


2k 


B 


(43) 


where  To  is  the  Doppler  cooling  limit.  For  Ca+,  the  Doppler  limit  is  0.53  mK. 

Key  to  Doppler  cooling  is  the  presence  of  a  fast  decaying  excited  electronic  state 
that  will  spontaneously  decay  to  the  ground  state.  In  addition,  it  is  convenient  if 
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the  transitions  are  accessible  with  commercially  available  lasers  [96,  100,  101].  To 
date  the  atomic  ions  that  have  been  cooled  include  Be+  [102],  Mg+  [103],  Ca+  [104], 
Ba+  [105],  Sr+  [106],  Yb+  [107],  Hg+  [108],  and  Cd+  [109].  With  the  exception  of  Be+ 
and  Mg+,  the  aforementioned  atomic  ions  do  not  have  a  simple  two-level  structure 
and  have  other  decay  channels.  These  are  are  plugged  with  secondary  re-pumping 
lasers  to  close  the  cycling  transition.  Ytterbium  possesses  several  transitions  that  are 
suitable  for  optical-frequency  reference  standards  (clocks)  such  as  the  2Si/2  O  2F7/2 
transition  at  467  mu. 

2.5  Laser  cooling  of  Yb+ 

2.5.1  Introduction 

The  ytterbium  ion  has  several  properties  that  make  it  an  ideal  candidate  for  many 
applications  in  physics  such  as  an  optical  reference  standard,  quantum  computing, 
and  investigations  of  time  variation  in  fundamental  physics  [110,  111],  Ytterbium 
ions  have  several  clock  transitions  such  as  the  ultra  narrow  2S'i/2  — >  2F\/2  transition 
at  467  nm,  the  2*Si/2  — >  2D5/2,  and  many  others  [107,  109,  112-115].  Ytterbium  also 
has  an  isotope  with  a  nuclear  spin  /  =  1/2  (iaYb)  with  high  abundance  (14.28%). 

2.5.2  Methods 

Ytterbium  ions  are  trapped  in  an  UHV  chamber  with  background  pressure  f»10~n 
torr  in  a  4.5  inch  octagon  chamber  pumped  by  a  50  L/sec  pump  (Daniway  pump).  The 
trap  is  5-segment  with  five  electrodes  and  a  radius  (r0  =  0.6 mm).  RF  voltage  (200V) 
driven  at  14.41  MHz  is  supplied  through  two  diagonalized  opposed  electrodes  while  the 
other  two  are  held  to  RF  ground.  DC  voltage  (~  25V)  is  supplied  on  the  axial  side.  An 
optical  seRip  showing  the  laser-layout  is  shown  in  the  schematic  below  in  Fig.  6.  Four 
lasers,  369  nm,  399  nm,  935  nm,  and  638  nm,  are  needed  to  trap  ytterbium  which 
will  be  explained  below.  About  20  mW  of  the  399  nm  laser  is  combined  with  about 
1  mW  369  nm  laser  with  a  bandpass  filter  (Semrock  395/20  nm)  that  reflects  the 
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399  nm  laser  and  transmits  the  369  nm  laser,  which  are  sent  along  the  trap  axis  and 
are  focused  with  a  150  mm  achromatic  lens  to  the  center  of  the  trap.  This  reduces 
the  spot  size  of  the  laser  beams  to  about  50  /irn.  Similarly,  about  10  mW  of  the 
935  nm  laser  is  combined  with  about  3  mW  of  the  638  nm  laser  via  a  bandpass  filter 
that  reflects  the  935  nm  laser  and  transmits  the  638  nm  laser.  The  combined  lasers 
are  also  focused  to  the  center  of  the  trap  with  a  150  mm  achromatic  lens  to  a  spot 
size  of  about  80  /mi. 


Filter  Trap  Dichroic  mirror 


Figure  6:  A  schematic  showing  the  optical  layout  of  the  lasers  to  trap  ytterbium 
ions.  The  lasers  are  coupled  into  the  trap  using  a  combination  of  filters  and  dichroic 
mirrors.  The  fluorescent  light  is  collected  with  a  setup  of  optics  and  focused  on  the 
EMCCD  camera  for  imaging  as  explained  in  Ref.  [116]. 
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Ytterbium  ions  are  produced  by  two  photon  ionization  of  neutral  ytterbium.  Neu¬ 
tral  ytterbium  is  introduced  in  the  chamber  by  resistively  heating  an  aluminum  oven 
filled  with  ytterbium  at  1.5  amps.  The  neutral  ytterbium  is  excited  on  the  1S'o  — >  lP\ 
transition  with  a  399  nm  laser  as  shown  in  Fig.  7a.  A  fluorescence  trace  of  neutral 
ytterbium  is  shown  in  Fig  7b.  A  second  photon  with  a  wavelength  shorter  than  394 
nm  is  enough  to  excite  the  electron  to  the  continuum.  Since  a  369  nm  laser  is  required 
for  Doppler  cooling,  it  is  also  used  for  the  second  photon  to  ionize  neutral  ytterbium. 


///////  / 
'  '  ' _  c  >  '  ' 

j  i 


394nm 


1 1 


399nm 


(a)  Yb  ionization  scheme  (b)  Yb  fluoresce 

Figure  7:  Ionization  of  neutral  ytterbium,  (a)  shows  the  isotopically  selective  ioniza¬ 
tion  scheme  from  neutral  Yb  to  the  continuum,  (b)  shows  the  fluorescence  of  neutral 
ytterbium  on  the  1S'o  — >  1P\  transition  imaged  at  399  nm. 

Doppler  cooling  of  Yb+  relies  on  the  cycling  transition  at  2Si/2  -H-  2P\/2  near 
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369  nm  as  shown  in  Fig.  8.  Occasionally  (~0.005%),  the  ion  will  decay  to  the  2D3/2 
meta  stable  state,  which  has  a  lifetime  of  52.7  ms  [117].  A  935  nm  laser  is  used  to 
pump  the  ion  from  the  2D3/2  — *  3  [3/2]  !/2 ,  which  spontaneously  decays  back  to  the 
2S 1/2  -H-  2 P\/2  cycling  transition.  Collisions  with  background  gas  causes  the  ions  to 
fall  in  to  the  long  lived  2Fj/2  state  which  has  a  lifetime  of  5.4  years  [118].  Another 
laser  at  638  nm  pumps  the  ion  to  the  1  [5/2]5/2  state,  which  decays  to  the  2D3/2  state 
and  then  back  to  the  cycling  transition.  The  odd  looking  spectroscopic  notation  of 
the  Yb+  states  are  caused  by  the  fact  that  those  states  are  no  longer  defined  by  the 
normal  L  —  S  coupling.  One  electron  from  the  core  /-shell  of  ytterbium  is  excited 
to  the  s  orbital  and  the  spin  of  the  s  are  coupled  with  angular  momentum  L.  While 
the  4/13  core  shell  remains  L  —  S  coupled  with  angular  momentum,  Jc.  The  sum  of 
the  two  momenta  K  —  L  +  Jc  ,  gives  the  resultant  coupling  which  is  shown  in  the 
square  brackets  [118]. 

2.5.3  Results 


Trapped  ytterbium  ions  are  shown  in  Fig.  8b.  The  ions  are  imaged  via  an  electron 
multiplied  charged  couple  device(EMCCD)  camera  (Princeton  Instruments  PI-532) 
with  a  lens  system  with  a  numeral  aperture  (N. A. =0.43)  and  a  magnification  of  10  x 
explained  in  detail  in  subsection  4.5  [116].  The  image  shows  a  Yb+  crystal  of  a  few 
hundred  ions.  In  order  for  the  ions  to  crystallize,  the  ratio  of  the  potential  energy  to 
the  kinetic  energy  must  be  greater  than  173  [119].  That  is, 


r 


1  q2 
47T6o  awksT 


>  173, 


(44) 


where  Y  is  the  Coulomb  coupling  parameter,  eo  is  the  permittivity  of  the  vacuum,  q  is 
the  charge  of  the  ion,  ks  is  the  Boltzmann  constant,  T  is  the  temperature,  and  aw  is 
the  Wigner-Seitz  radius  defined  as  47ra3 /3  =  1  /n0  with  n0  being  the  particle  density. 
The  Wigner-Seitz  radius  describes  the  effective  distance  between  the  ions. 
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(a)  Yb+  energy  levels  (b)  A  crystal  of  trapped  Yb+  ions 


Figure  8:  Energy  levels  for  laser  cooling  a  Yb+  ion  (a).  The  2S\/2  -H-  2P\/2  transition 
is  used  for  Doppler  cooling.  However,  Yb+  has  two  meta-stable  states  that  the  ion 
decays  to  and  need  to  be  re-pumped;  one  at  935  nm  and  another  at  638  nm.  A  crystal 
of  trapped  Yb+  ions  is  shown  in  (b). 

We  were  able  to  trap  different  isotopes  of  Yb+  (Table  1).  We  identified  the  differ¬ 
ent  isotopes  comparing  our  frequency  with  the  expected  frequency  values  for  those 
transitions  (Ref.  [113,  117,  120]).  Due  to  hyperfine  splitting,  the  171  isotope  requires 
are  microwave  source  of  12.643  GHz  to  depopulate  the  ions  from  the  state  that  is  not 
in  the  cycling  level.  To  demonstrate  that  we  had  trapped  171Yb+,  we  turned  on  and 
off  the  microwave  source  and  the  ions  appeared  and  disappeared  accordingly. 
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Table  1:  Trapping  frequencies  for  different  Yb+  isotopes.  The  percentages  in  brackets 
reflect  the  natural  abundance  of  each  isotope.  The  frequencies  were  measured  by 
the  lab  wavemeter  (HighFinesse  WS7).  The  frequencies  are  reported  in  THz.  The 
right  most  column  shows  corresponding  frequencies  for  the  neutral  Yb  atoms  for  the 
various  isotopes.  Other  isotopes  of  Yb,  were  not  observed  as  168  has  very  low  abun¬ 
dances  and  173  has  many  more  decay  channels  that  required  additional  re-pumping 
lasers.  The  2P7/ 2  -H-  2[5/2]!/2  transition  was  measured  to  be  469.43921  THz  and  no 
isotopic  difference  was  observed.  This  was  verified  by  blocking  the  638  nm  light  and 
waiting  for  the  ion  to  go  dark.  The  laser  would  then  be  unblocked  and  the  ion  would 
appear  again. 


Yb+  (amu) 

2Sl/2  O  2P1/2(THz) 

2D3, 2  ++  2[3/2]1/2(THz) 

Yb  -H-  1  Pi  (THz) 

176  (12.76%) 

811.2902 

320.5746 

751.52706 

174  (31.83%) 

811.2915 

320.5720 

751.52674 

172  (21.83%) 

811.2922-30 

320.56940 

751.52650 

170  (3.04%) 

811.2939-44 

320.5657 

751.52612 

171  (14.28%) 

811.28854 

320.568898 

751.52610 

2.6  Laser  cooling  Ca+ 

The  rest  of  the  experiments  described  in  this  thesis  use  40Ca+  for  laser  cooling.  40Ca+ 
ions  are  produced  from  two-photon  ionization  of  neutral  calcium  atoms.  A  cloud  of 
calcium  atoms  is  produced  by  resistively  heating  an  oven  filled  with  calcium.  The 
calcium  is  the  ionized  first  by  exciting  the  1S'o  1P\  transition  with  a  423  nm 
laser.  Then  a  second  photon  around  388  nm  ionizes  the  neutral  calcium  as  shown 
in  Fig.  9a.  Scanning  the  423  nm  laser  shows  a  fluorescence  of  the  Doppler  shifted 
excitation  of  the  neutral  calcium  atoms  as  shown  in  Fig.  9b. 
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(a)  Ca  ionization 
scheme 

Figure  9:  The  ionization  scheme  for  neutral  calcium,  (a)  Isotope  selective  ionization 
of  calcium  is  performed  by  exciting  the  neutral  calcium  along  the  1So  — >  lP\  transition 
with  a  tuned  423  nm  photon  and  another  photon  with  wavelength  of  about  388  nm 
ionizes  it  to  the  continuum,  (b)  An  image  of  fluorescent  neutral  calcium  atoms  imaged 
on  the  EMCCD  camera. 

The  Ca+  ion  has  a  ground  state  of  4  2Si/2-  Two  hyperhne  levels  lie  above  the 
ground  state,  3  2D3/2  and  3  2D5/2.  Above  these  first  excited  states  are  another  two 
hyperhne  levels,  4  2Pi/2  and  the  4  2P3/2  as  shown  in  Fig.  10a.  Doppler  Cooling  is 
realized  by  cycling  on  the  397  nm  2S'i/2  2Pi/2  transition  as  shown  in  Fig.  10a, 

which  has  a  lifetime  7.1  ns.  Occasionally,  one  in  twelve  cycles,  the  ion  decays  to  the 
D3/2  electronic  level,  which  has  a  lifetime  of  1  s,  and  is  re-pumped  with  a  866  nm 
laser  along  the  2D3/2  2P3/2  transition  back  to  the  cycling  transition.  To  image  the 

trapped  ions,  fluorescence  from  the  397  nm  transition  is  collected  with  an  objective 


50  pm 


(b)  Ca  fluorescence 


27 


lens  (magnification=10x)  as  shown  below  in  Fig.  10b. 


(b)  A  chain  of  40Ca+  ions 


Figure  10:  Energy  levels  for  laser  cooling  the  40Ca+  ion  (a).  The  2Si/2  2P\/2 

transition  is  used  for  laser  cooling  near  397  nm.  An  additional  laser  at  866  nm  laser 
is  used  to  re-pump  the  ion  from  2-D3/2  level.  The  729  nm  and  854  nm  lasers  are  used 
for  sideband  cooling,  which  is  explained  in  detail  in  section  2.7.  A  chain  of  cooled 
40Ca+  ions  (b)  imaged  on  an  EMCCD  camera  is  shown. 


The  frequency  of  the  transitions  is  measured  using  a  calibrated  commercial  Fabry- 
Perot  Interferometer  called  the  WS7  (HighFinesse).  The  accuracy  of  the  measured 
transition  frequencies  is  limited  by  the  WS7,  which  has  an  absolute  accuracy  of  60 
MHz.  The  measured  transitions  are  shown  in  the  table  below  (Tab.  2). 


Table  2:  Frequencies  for  different  transitions  of  40Ca+  measured  in  the  laboratory 
with  the  wavemeter  (HighFinesse  WS7).  The  frequencies  are  those  observed  in  air 
and  the  right  most  column  is  the  frequency  of  the  neutral  calcium  transition. 


40Ca+  Transition 

Frequency  (THz) 

40 Ca  1S'0  O  1P1(THz) 

2Si/2  -H-  2 P\/2 

755.2227 

709.0783 

2S\/2  2  71,5/2 

411.0421 

2_D3/2  2 P\/2 

346.0003 

2P>5/2  ^  2 P3/2 

352.6734 

2. 7  Sideband  cooling 

Narrow  transitions  on  sidebands  of  trapped  atomic  ions  are  ideally  suited  for  precision 
spectroscopy  of  molecular  ions  [121].  Before  precision  spectroscopy  can  be  achieved, 
the  ion-molecular-ion  pair  need  to  be  sideband  cooled  to  remove  the  motional  energy 
of  the  ion  in  the  trap.  For  sideband  cooling,  the  trap  frequency  must  be  larger  than 
the  laser  linewidth  and  the  decay  rate  of  the  excited  state.  This  makes  it  possible 
to  address  each  individual  sideband,  which  is  critical  for  sideband  cooling.  When  the 
laser  is  detuned  to  the  red  sideband,  the  sideband  with  lower  motional  energy,  an 
excitation  to  the  upper  level  results  in  a  reduction  of  motional  energy. 

Since  the  sideband  cooling  state  is  long-lived,  it  is  coupled  to  a  short-lived  state 
with  another  laser  for  faster  cooling  rate.  Spontaneous  decay  does  not  change  the 
vibrational  quantum  number  of  the  ions.  For  Ca+  the  2Si/2  O  2D5/2  transition  is 
used  for  sideband  cooling  with  a  729  nm  laser.  The  long  lived  2D5/2  state  is  coupled 
to  the  short-lived  2P3/2  state  with  an  854  nm  laser  to  increase  the  cooling  rate.  A 
magnetic  held  splits  the  degenerate  level  in  the  sideband  transition  to  about  ten 
states.  For  sideband  cooling,  the  2S\/2{m  =  —1/2)  — »  2D5/2(;m  =  —5/2)  since  there 
is  no  Zeeman  level  lower  that  the  2D^/2((m  =  —5/2)  state.  Hence,  there  will  not  be 
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a  sideband  closer  to  energy  with  any  of  the  red  sidebands  of  the  2-D5/2(m  =  —5/2) 
state.  It  is  convenient  to  start  the  theoretical  description  of  sideband  cooling  from  the 
Lamb-Dicke  regime  since  that  is  where  sideband  cooling  experiments  begin;  that  is, 
sideband  cooling  is  proceeded  by  Doppler  cooling.  The  Lamb-Dicke  regime  is  when  the 
ions  are  confined  in  a  space  much  smaller  than  the  wavelength  of  the  transition.  The 
Lamb-Dicke  regime  is  characterized  by  the  Lamb-Dicke  parameter  rj  which  is  defined 
as, 

rj  =  kx o  <<  1,  (45) 


where  k  is  the  wavevector  and  Xq  =  yjh/2mu  and  u  describes  the  harmonic  oscillation 
frequency.  Furthermore, 


2  (hk  cos  6)2  . 


(46) 


where  6  is  the  angle  between  the  wave  vector  and  the  trap  axis.  Thus,  the  square  of 
the  Lamb-Dicke  describes  the  ratio  between  photon  recoil  energy  and  the  energy  level 
spacing  the  in  the  harmonic  oscillator.  In  the  Lamb-Dicke  regime,  the  interaction 
Hamiltonian  can  be  written  as, 


(0i,  n' |  Hi  | -0o,™) 


^he~iAtn(n' \ei71^+a) 


n), 


(47) 


where  n  =  0,1,2....  are  the  levels  of  the  harmonic  oscillator  and  Q  is  the  Rabi  fre¬ 
quency.  Because  rj  is  very  small,  the  exponential  term  in  the  Hamiltonian  can  be 
expanded  using  the  Taylor  series  as, 


ei’<a,+“>  =  1  +  ir)(a<  +  a)  +  0(i;2). 


(48) 


This  results  in  processes  that  change  the  vibrational  quantum  number  by  more  than 
one  to  be  suppressed  and  favor  changes  by  ±1.  The  interaction  Hamiltonian  defines 
the  allowed  transitions  which  broadly  fall  into  three  categories: 

•  Red  sideband  transition- These  are  transitions  where  the  change  in  vibrational 
quantum  number  is  An  =  —1.  That  is  the  transition  results  in  an  overall 
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decrease  in  the  motional  energy. 


^n— i,n  =  ??\MV  (49) 

•  Carrier  transition- There  are  transitions  where  excitation  to  a  higher  electronic 
state  results  in  no  change  in  the  vibrational  quantum  number  i.e.  An  =  0  and, 

^n,n  =  ^o(l  ^rfn).  (50) 

•  Blue  sideband  transition-  These  are  transitions  where  an  electronic  change  in¬ 
creases  the  vibrational  quantum  number  by  An  =  +1  and, 

&n+i,n  =  V'/n  +  l^o-  (51) 

The  Rabi  frequencies  of  the  sidebands  are  proportional  to  rf.  The  transitions  are 
shown  schematically  in  Fig.  11. 
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Figure  11:  An  illustration  of  the  driven  transition,  (a)  shows  the  red  sideband  transi¬ 
tion,  (b)  shows  the  carrier  transition,  and  (11c)  shows  the  blue  sideband  transition,  g 
and  e  represent  the  electronic  state,  while  0,  1  and  2  represent  the  lowest  vibrational 
states. 


Repeated  transitions  on  the  red  sideband  removes  the  vibrational  quanta  to  the 
motional  ground  state.  This  process  is  called  sideband  cooling  and  is  used  to  cool  ions 
beyond  the  Doppler  cooling  limit  [122-130]. 

2.8  Non- destructive  identification  of  molecular  ions 

To  do  chemistry  in  ion  traps,  the  identity  of  the  reactant  must  be  ascertained  accu¬ 
rately.  Because  of  the  controlled  nature  of  chemical  reactions  in  ion  traps,  only  mass 
spectrometry  is  sufficient  to  accurately  identify  ions  inside  the  trap.  Traditionally, 
mass  spectrometry  was  carried  out  with  time  of  flight  (TOF)  and  quadrupole  mass 
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filter  (QMF)  techniques.  These  methods  destroy  the  ions  on  measurement.  An  alter¬ 
native  technique  would  be  to  measure  the  mass  of  the  ion  in  situ  without  destroying 
it  so  that  experiments  can  be  carried  out  on  that  ion  after  identification.  We  used  a 
non-destructive  technique  that  relates  the  oscillation  frequency  of  the  ions  in  the  trap 
with  their  mass  [131,  132],  Since  the  ions  are  coupled  by  their  Coulombic  interaction, 
any  perturbation  on  one  ion  affects  the  other  ion.  Thus  for  this  technique  to  work 
the  identity  and  frequency  of  one  ion  must  be  known  prior  the  addition  of  the  second 
unknown  ion.  For  two  ions,  the  axial  frequencies  are  governed  by  the  equation 

^+/_  =  [(1  +  fj)±  y/1  -h-  +  h-2K2,  (52) 

where  u+  is  the  breathing  mode  (out-of-phase  oscillation)  and  v_  is  the  center  of  mass 
mode  (in-phase  oscillation)  of  the  oscillation  of  the  ions,  /i  is  the  ratio  of  the  two  ions 
Mi/ M2.  To  identify  the  resonance  frequency  of  the  ions,  a  sinusoidal  varying  voltage 
is  applied  on  one  of  the  electrodes.  When  the  frequency  of  the  applied  voltage  is  on 
resonance  with  the  secular  frequency  of  the  ion  crystals,  the  crystal  mode  motion  is 
excited  and  the  fluorescence  signal  of  the  lase-cooled  atomic  ions  is  perturbed.  This 
will  show  up  as  either  a  clip  or  a  spike  in  the  fluorescence  spectrum  depending  on  the 
detuning  of  laser.1  Monitoring  the  fluorescence  signal  as  a  function  of  the  frequency 
of  the  drive  voltage,  the  identity  of  the  unknown  ion  can  be  deduced  from  the  above 
equation.  This  technique  was  pioneered  by  the  Baba  and  Waki  group  in  Japan  and 
the  Drewsen  group  in  Denmark  [61,  131-133].  An  example  using  this  technique  is 
shown  below  in  Fig.  12.  This  method  is  commonly  referred  to  as  a  tickle  scan  because 
the  voltage  tickles  the  secular  motion  of  the  ions  in  the  trap. 

1If  the  laser  is  red-detuned  from  the  atomic  transition,  an  excitation  of  the  oscillation  mode  by 
the  applied  voltage  causes  a  heating  of  the  crystal,  which  brings  the  red-detuned  laser  on  resonance 
with  the  laser-cooled  ion,  resulting  in  an  increase  in  fluorescence  counts,  a  peak.  However,  if  the 
laser  is  on  resonance  with  the  laser-cooled  atomic  ion  prior  to  the  excitation,  the  heating  caused  by 
the  excitation  of  the  oscillation  mode  will  result  in  overall  heating  of  the  trapped  crystal,  leading  to 
poor  cooling  and  low  fluoresce  counts.  This  shows  up  as  a  dip  in  the  spectrum. 
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Figure  12:  The  secular  frequency  scan  on  the  center  of  mass  mode  showing  two  peaks. 
The  circles  are  the  data,  while  the  lines  are  Gaussian  fits  to  the  data.  The  peak  at 
359  kHz  (red  curve)  shows  a  tickle  scan  of  one  40Ca+  ion  where  else  the  peak  at  348 
kHz  (blue  curve)  shows  a  tickle  scan  of  one  40Ca+  and  another  unknown  ion.  The 
identity  of  the  unknown  ion  was  determined  to  be  44Ca+  based  on  equation  52. 


In  this  example,  two  species  of  Ca+  isotopes  were  trapped.  One  was  the  known 
40Ca+,  while  the  other  was  unknown.  Prior  to  trapping  two  species,  only  one  40Ca+ 
was  trapped.  The  secular  frequency  of  the  known  40Ca+  was  determined  by  applying 
an  oscillating  frequency  on  the  electrodes  while  the  laser  was  slightly  red  detuned 
from  the  cooling  transition.  An  excitation  of  the  secular  frequency  of  the  trapped  ion 
results  in  heating  of  the  ion  which  causes  the  Doppler  shifted  laser  to  be  closer  to 
resonance  of  the  atomic  transition  increasing  the  cooling.  This  shows  up  in  the  spectra 
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as  an  increase  in  the  fluorescence  counts.  The  red  trace  in  Fig.  12  is  a  Gaussian  fit  to 
the  data  that  determined  the  frequency.  Another  isotope  of  Ca+  was  trapped  along 
with  the  40Ca+  ion.  However,  this  isotope  was  not  fluorescent,  hence  its  identity  was 
unknown.  A  tickle  voltage  was  applied  on  the  crystal  and  the  resultant  scan  is  shown 
as  the  solid  circles.  The  data  was  fit  with  the  blue  Gaussian  line.  Using  equation  52, 
the  identity  of  the  unknown  isotope  was  found  to  be  44Ca+.  The  determined  mass 
was  within  0.5  amu  difference  with  the  predicted  mass. 
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CHAPTER  III 


MOLECULAR  SPECTROSCOPY  AND  40CaH+ 
MOLECULAR  ION 


To  assign  the  spectra  of  a  molecule,  several  properties  of  the  molecule  have  to  be 
known  first.  These  properties,  called  spectroscopic  constants,  come  from  the  Schrodinger 
equation  describing  the  molecular  system.  In  this  chapter,  the  basics  of  spectroscopic 
transitions  will  be  reviewed  as  well  as  a  summary  of  the  40CaH+  molecular  ion. 

3.1  Spectroscopy  and  selection  rules 

The  molecular  Hamiltonian  has  many  degrees  of  freedom  and  can  be  written  as  [134, 
135], 


where  Mk  is  the  mass  of  the  nucleus,  me  is  the  electron  mass,  £o  is  the  permittivity 
of  free  space  constant,  Z%  is  the  charge  of  the  nucleus,  e  is  the  charge  of  the  electrons, 
R  is  the  distance  between  the  two  nuclei  and  rt  is  the  coordinate  of  the  electrons. 
The  first  term  describes  the  kinetic  energy  of  the  nuclei,  the  second  term  the  kinetic 
energy  of  the  electrons.  The  terms  in  the  square  brackets  describe  the  potential  energy 
generated  by  the  repulsion  of  the  nuclei,  electron-electron  repulsion  and  electron- 
nuclear  attraction  respectively.  Because  the  motion  of  the  nuclei  is  much  slower  than 
that  of  the  electrons,  the  nuclei  appear  static  compared  to  the  electronic  motion.  The 
Born-Oppenhcimer  approximation  [136]  allows  for  the  separation  of  nuclear  motion 
and  electronic  motion.  Thus  the  wavefunction  can  be  separated  as, 


1 molecule')  nuclear)  %  | ^electronic)- 
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Thus  the  Hamiltonian  can  be  separated  into  the  electronic  contribution,  the  vibration 
and  rotation  of  the  nuclei  as  follows: 


H, 


molecule 


=  H, 


electronic 


H, 


vibrational 


H, 


rotational  • 


(55) 


The  Hamiltonian  of  the  vibrational  and  rotational  contributions  for  a  diatomic  molecule 
can  be  approximated  to  first  order  as, 


1 

(56) 

H vibrational  ^ e\y  “1“ 

H rotational  B  J (J  H-  1), 

(57) 

where  the  constants  expressed  in  (cm  x)  are, 


1  *  k 

We  = 

IlC  y  yU 


B  =  h1 
he' 21' 


(58) 

(59) 


where,  uje  is  the  angular  frequency  of  the  transition,  c  is  the  speed  of  light,  k  is 
the  force  constant,  h  is  Planck’s  constant,  h  is  the  reduced  plank  constant,  /  is  the 
molecule’s  inertia  in  terms  of  the  reduced  mass, 


i  =  /*>'). 


(60) 


mirri2 

mi  +  m‘2  ’ 


(61) 


with  r0  being  the  equilibrium  bond  length.  Since  vibrational  and  rotational  motions 
are  coupled,  the  energy  of  the  molecule  is  expressed  as  the  sum  of  the  rovibrational 
energy, 


f^rovib  d  "f  1). 


(62) 


The  above  equations  describe  ideal  situations  of  a  perfect  harmonic  oscillator  and 
rigid  rotor.  In  reality,  the  energy  of  a  molecule  becomes  perturbed  at  high  energies. 
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Thus  corrections  are  added  to  account  for  the  anharmonicity  of  a  molecule.  An  ap¬ 
proximation  to  this  anharmonicity  is  the  Morse  potential  and  it  takes  the  form, 

V{r)  =  De{(  (63) 

where  De  is  the  depth  of  the  potential  i.e.  the  dissociation  energy  for  the  molecule 
and  a  =  \Jke/2De.  The  energies  of  the  rovibrational  oscillator  become  second  order 
in  v  and  fourth  order  in  J  as, 

X  1 

Erovib  =  +  — )  —  uexe(y  +  — )2  +  BVJ(J  +  1)  —  DVJ\J  +  l)2.  (64) 

The  new  terms  ujexe  and  D  are  the  anharmonicity  constant  and  the  centrifugal  distor¬ 
tion  constant  respectively.  In  order  to  have  a  rovibrationl  spectrum,  there  has  to  be 
a  change  in  the  dipole  moment  accompanying  that  vibration  or  rotation.  The  change 
in  vibrational  motion  for  a  harmonic  oscillator  is  governed  by  the  selection  rule  that 
A v  =  ±1.  However,  due  to  the  anharmonicity  of  molecular  vibrations  other  tran¬ 
sitions  such  as  A v  =  ±2,  ±3,  ±4...  are  possible  albeit  at  decreasing  intensities.  On 
the  other  hand,  for  molecular  rotations,  the  selection  rule  A  J  =  0,  ±1  holds  true 
even  when  rotor  is  distorted.  This  is  because  D  <  B  hence  the  distortion  is  very 
small.  Transitions  with  A  J  =  0  are  observed  when  there  is  a  change  in  electronic  en¬ 
ergy.  Transitions  with  A  J  =  +1  are  labeled  the  R  branch,  transitions  with  A  J  =  0  are 
called  the  Q  branch,  and  transitions  with  A  J  =  —  1  are  called  the  P  branch.  Thus 
from  an  observed  rovibrational  spectrum,  we  can  calculate  the  vibrational  and  ro¬ 
tational  constants,  oje,  ooexe,  Bv,  and  Dv.  Conversely,  if  we  can  calculate  from  first 
principles  the  rovibrational  constants,  we  can  approximately  predict  where  the  lines 
for  that  spectrum  will  be. 

When  spectra  is  accompanied  with  changes  in  electronic  energy,  electronic  transi¬ 
tion  selection  rules  apply  together  with  the  rovibrational  selection  rules.  These  rules 
are  governed  by  the  angular  momentum  of  the  electrons  about  the  internuclear  axis. 
The  motion  of  the  electrons  about  the  axis  gives  rise  to  orbital  angular  momentum 
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L  which  precesses  about  the  nuclear  axis  and  its  projection  about  the  axis  is  called 
Ml.  At  faster  precession  speeds.  L  ceases  to  be  a  good  quantum  number,  while 
its  projection  remains  a  good  quantum  number.  Hence  electronic  states  are  clas¬ 
sified  based  on  Mg,  which  is  labeled  A  =  M L \ .  The  values  that  A  can  take  are 

A  =  0, 1, 2,  3 . L.  The  molecular  states  are  labeled  accordingly  to  each  value  of  A  as 

E,  n,  A,  <f>  and  so  on.  The  molecules  with  E  state  has  symmetry  properties  about  the 
plane  through  the  internuclear  axis  and  can  take  the  form  E+  or  E“.  The  electrons 
in  the  molecules  have  spin  S  ,  which  can  also  be  projected  on  to  the  internuclear  axis 
where  it  is  called  Mg.  Similar  to  the  case  orbital  angular  momentum,  S  no  longer 
becomes  a  good  quantum  number  while  its  projection  remains  defined.  For  bookkeep¬ 
ing,  we  call  the  projection  Mg  to  be  E  in  molecules.  There  are  2S  +  1  possible  values 
of  E.  The  total  angular  momentum  of  electronic  motion  becomes  the  sum  of  the 
orbital  angular  momentum  and  the  spin  angular  momentum; 

H=|A  +  E|,  (65) 

Where  can  be  |A  —  S  | ,  |  A  —  E|-t-l,...|A  +  E|. 

When  electronic  change  is  coupled  with  rotational  change,  the  two  motions  become 
coupled  which  give  rise  to  changes  in  the  spectra.  There  are  different  ways  in  which 
the  momenta  can  couple  and  the  different  schemes  are  grouped  into  Hund’s  couple 
cases.  While  the  coupling  cases  range  from  (a)  to  (e),  the  most  common  cases  are 
Hund’s  coupling  cases  (a)  and  (b),  which  will  be  reviewed  here.  For  the  rest  the 
reader  can  refer  to  [1]. 
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3.2  Hund’s  case  (a) 


Figure  13:  A  vector  diagram  showing  Hund’s  case  (a).  It  is  defined  by  strong  spin- 
orbit  coupling  to  the  internuclear  axis  to  give  the  resultant  angular  momentum 


In  Hund’s  case  (a),  the  motion  of  electrons  couples  strongly  to  the  internuclear  axis 
to  give  a  large  spin-orbit  coupling.  N  represents  the  rotation  of  the  nucleus.  N  and 
precess  around  J,  which  is  the  resultant  angular  momentum.  The  rotational  energy 
of  the  molecule  becomes, 

E(J)  =  B[J(J +1) -n2}  (66) 
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3.3  Hund’s  case  (b) 


N 


Figure  14:  A  vector  diagram  showing  Hund’s  case  (b).  The  spin  angular  momentum 
is  weakly  coupled  to  the  internuclear  axis,  instead  it  couples  to  the  sum  of  the  nuclear 
rotation  and  electronic  angular  momenta 

When  A  =  0  and  S  /  0,  the  spin  is  not  strongly  coupled  to  the  internuclear  axis. 
Under  this  condition,  is  no  longer  a  good  quantum  number.  Instead  A  and  N  form 
the  resultant  angular  momentum  K  without  the  spin.  The  spin,  S  couples  with  K 
to  form  the  resultant  angular  momentum  J.  For  CaH+  the  ground  state  is  described 
by  rigid  rotor  since  A  =  0  and  S=  0,  some  of  the  excited  electronic  states  by  Hund’s 
case  (b). 

3.4  Electronic  selection  rules 

For  molecules  where  A  is  well  defined,  the  selection  rule  is  that  transitions  must  have 
A  =  0,  ±1.  Thus  the  possible  electronic  transitions  are  E  —  E,  E  —  n  and  so  on.  I11 
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addition,  E  states  cannot  change  symmetry  during  a  transition.  So  the  transitions 
E+  yy  E+  and  E“  yy  E~  are  possible  but  not  E+  yy  E-. 

3.5  40  CaH+ 

Similar  to  other  alkaline  earth  hydrides  ions,  40CaH+  molecular  ions  are  of  particular 
scientific  interest  because  of  potential  applications  to  astrophysics  and  fundamental 
physics.  Although  not  detected,  these  hydrides  are  thought  to  be  present  in  ISM 
and  comets  [14,  15,  20,  137-140].  Kajita  and  Moriwaki  have  identified  halides  of 
alkali  metals  as  ideal  candidates  for  the  testing  of  //,  proton-to-electron  mass  ratio. 
40CaH+  has  been  proposed  as  a  good  candidate  for  these  measurements  [141],  Kajita 
and  coworkers  have  also  performed  high  precision  ab  initio  calculations  on  40CaH+ 
vibrational  overtones  [25].  In  particular,  we  are  interested  in  using  our  laboratory 
techniques  to  perform  the  first  measurements  of  40CaH+  vibrational  spectra.  We  will 
focus  on  the  transitions  from  v'  —  9  4—  v  —  0  and  v'  —  10  <—  v  —  0.  These  occur  at 
883.9  nm  and  813.8  nm,  respectively. 

Despite  the  lack  of  experimental  evidence,  there  has  been  a  significant  amount  of 
theoretical  investigations  into  the  40CaH+  molecule.  To  measure  the  spectra  of  an 
unknown  molecule,  we  consulted  the  theory  calculations  as  a  guide.  Table  3  shows  a 
summary  of  the  spectroscopic  constants  found  in  the  literature. 
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Table  3:  A  summary  of  published  theoretical  spectroscopic  constants  of  40CaH+  found 
in  literature.  For  the  experiments  described  in  this  thesis,  we  use  the  values  from  [24, 
25].  The  other  calculations  used  a  lower  level  of  theory,  while  the  others  did  not 
evaluate  the  vibrational  wavefunctions. 


Reference 

D0  (cm  4) 

ae(cm  4) 

r(A) 

Be  (cm  4) 

tce(cm  4)  ccexe(cm 

Ref 

[17] 

14809 

0.13 

1.936 

4.609 

1511.0 

23.5 

Ref 

[14] 

2.085 

1416.8 

Ref 

[20] 

17000 

1.864 

4.85 

1504 

21.01 

Ref 

[139] 

16719.9 

0.0985 

1.881 

4.85 

1482 

21.01 

Ref 

[142] 

15500 

1.926 

1482 

Ref 

[21] 

17972 

1.873 

1453 

Ref 

[22] 

17876 

1.86 

Ref 

[143] 

17876.80 

1.936 

4.609 

1468 

Ref 

[24] 

14804 

1.936 

4.609 

1511 

Ref 

[25] 

16267.39 

0.09 

1.926 

4.611 

1457.7 

The  energetics  of  the  40Ca+  +  H2  reaction  are  similar  to  those  of  neutral  Ca  and 
H2  [55,  139,  144-147]. 1  The  bond  energy  of  H2  is  4.478  eV,  while  the  bond  energy  of 
40CaH+  is  1.788  eV.  This  makes  the  reaction  endogernic  with  an  energy  difference  of 
2.77  eV  [144],  To  get  over  the  barrier,  the  40Ca+  is  excited  with  the  397  nm  laser  (3.12 
eV)  to  the  2Pi/2  state  in  order  for  the  reaction  40Ca+  +  H2— »  40CaH+  +  H  reaction 
to  occur.  Once  reacted,  the  molecule  forms  a  stable  compound  with  an  equilibrium 
bond  length  of  about  1.92  A.  The  molecule  forms  an  ionic,  closed  shell  bond  of  Ca2+- 
H~  with  a  1/R  curve.  The  ground  state  of  the  molecule  has  an  asymptote  40Ca+(4s) 

lrThis  is  because  the  energy  barrier  for  the  neutral  Ca  reaction  with  H2  is  close  in  energy  to  that 
of  ionized  calcium  with  hydrogen.  Furthermore,  studied  reactions  between  Mg+  and  H2  had  similar 
reaction  dynamics  with  those  of  the  neutral  magnesium  and  hydrogen  [148-151]. 


43 


+  H(ls)  with  a  (4s  a)2  1£+  symmetry.  The  potential  energy  surfaces  of  40CaH+ 
are  shown  in  Fig.  15.  The  first  excited  state  has  an  asymptote  of  40Ca+(3d)  +  H(ls) 
with  the  electronic  configuration  of  4sa3dcr,  which  produces  the  molecular  states  1,3£, 
1,3II,  and  1,3A.  The  are  other  higher  electronic  state  which  are  below  the  ionic  limit, 
however,  for  this  thesis  only  the  first  excited  asymptote  will  be  considered. 


Figure  15:  Potential  energy  curves  of  the  ground  state  and  the  lowest  excited  states 
of  40CaH+  computed  at  the  EOM-CCSD/cc-pCVQZ  level  of  theory  provided  by  [152], 


The  40CaH+  molecule  possesses  a  permanent  dipole  moment,  which  makes  it  possi¬ 
ble  have  a  pure  rotational  spectrum  for  future  precision  measurements.  The  computed 
dipole  moment  is  shown  in  Fig.  16. 
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Figure  16:  The  dipole  moment  of  the  ground  state  of  10CaH+  provided  by  [141]. 

For  experimental  convenience,  we  chose  to  investigate  the  i/  =  9  v  =  0  and  the 
v'  =  10  v  =  0  expected  around  883  nm  and  813  nm  respectively.  This  is  because 
those  transitions  are  in  the  range  of  a  commercially  available  titanium-sapphire  laser 
without  any  modification  as  explained  in  detail  in  subsection  4.2.  Because  these  are 
higher-order  overtones  the  coupling  strength  becomes  very  weak  and  the  intensity 
becomes  very  weak.  The  expected  transition  strength  were  calculated  with  the  LEVEL 
program  to  compute  the  expected  matrix  elements  of  the  transitions  as  shown  in 
table  4  [153]. 
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Table  4:  Matrix  elements  of  the  Xu  •(—  Xu'R( 0)  transitions  calculated  with  the 
LEVEL  program  provided  by  Ref.  [154],  Here  the  v  indicates  the  lower  level  while 
u'  indicates  the  upper  level.  E(2)-E(l)  is  the  difference  in  energy  between  the  two 
states.  The  calculations  were  based  on  the  40CaH+  potential  given  by  [25]. 


dJ(J”) 

u-u' 

E(lower) 

E(2)-E(l) 

A  (Einstein) 

F-C  Factor 

(uJ\M\u'J') 

R(0) 

9-0 

-16714.06 

11320.89 

2.14694E-4 

1.21794E-16 

-3.7622E-5 

R(o) 

10-0 

-16714.06 

12295.00 

5.54952E-5 

2.28182E-17 

1.69003E-5 

3.6  Multiphoton  techniques  in  spectroscopy 

Because  many  spectroscopic  transitions  of  interest  generally  do  not  exactly  fall  in  the 
range  of  commercially  available  lasers,  it  is  a  common  tool  in  spectroscopy  to  use  a 
combination  of  laser  to  access  the  transition(s)  of  interest.  Multiphoton  techniques 
employ  the  absorption  of  two  or  more  photons  of  different  or  same  frequency  to  study 
transitions  that  may  not  be  easily  accessible  with  one  laser  or  the  detection  of  that 
state  is  is  not  accessible  with  other  detection  methods  such  as  fluoresce  quenching. 
Typically,  one  photon  excites  the  molecule  to  an  intermediate  state  and  then  a  second 
photon  excites  it  to  the  final  state.  From  these  techniques  several  properties  of  the 
molecule  can  be  studied.  If  the  first  photon  is  resonant  with  a  fluorescent  low  excited 
state,  the  second  transition  can  be  probed  by  scanning  the  second  laser  light  and 
observing  the  change  in  signal.  In  principle  many  lasers  can  be  employed  covering  a 
big  energy  gap. 

For  neutral  molecules,  typical  methods  of  detection  include  observing  the  change  in 
fluorescence  when  stimulated  emission  occurs  from  the  intermediate  state  once  second 
laser  is  turned  on,  ionizing  the  molecule  and  counting  the  number  of  ions  produced 
at  a  given  wavelength,  (e.g.  resonance  enhance  multiphoton  ionization(REMPI)). 
This  method  is  advantageous  because  the  observed  transitions  are  those  that  are 
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are  in  the  intermediate  state.  Therefore,  even  in  molecules  with  many  electronic 
states  with  overlapping  vibrational  states  can  be  discriminated  with  multiphoton 
spectroscopy.  Multiphoton  spectroscopy  has  been  applied  in  the  study  of  Doppler 
free  spectroscopy  in  sodium  [155],  hydrogen  bonding  in  aromatic-X  clusters  [156], 
rotational  spectroscopy  in  CaH  and  CaD  [157],  spectroscopy  of  biomolecules  such  as 
tryptophan  [158],  torsional  barriers  in  nonrigid  molecules  [159],  studies  of  Rydberg 
states  [160],  the  development  of  laser  photoelectron  spectroscopy  [161]  and  many 
other  experiments. 

3.6.1  Multiple  photon  dissociation 

Molecular  ions,  on  the  other  hand,  present  a  challenge  to  some  multiphoton  techniques 
since  electrostatic  repulsions  between  ions  limits  densities.  In  addition  the  methods  of 
detection  for  neutrals  become  difficult  for  ions.  In  REMPI,  ionization  is  used  to  probe 
the  intermediate  state  [162-164],  Molecular  ions  have  very  high  ionization  potentials 
that  would  be  impractical  for  the  most  part.  Most  molecular  ions  do  not  fluoresce  and 
makes  it  difficult  to  detect  them.  For  these  reasons,  resonance  enhanced  multiphoton 
dissociation  (REMPD)  is  employed  in  the  study  of  molecular  ions  [81,  165-169].  A 
schematic  of  REMPD  is  shown  in  figure  17. 
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Dissociation 


the  dissociative  state  Dx  via  an  intermediate  state  Si  from  the  ground  state  S0  by 
absorbing  two  photons  with  energy  Eq  and  E\ .  The  spontaneous  emission  rate  Kr, 
competes  with  the  REMPD  process  but  because  the  REMPD  process  is  significantly 
faster  than  spontaneous  emission,  it  is  neglected  in  the  theoretical  treatment. 

In  REMPD,  a  molecular  ion  in  the  ground  state  is  excited  to  an  immediate 
state  (rovibrational,  electronic)  by  one  laser  v\.  A  second  laser,  z/2,  couples  that 
excited  state  with  dissociative  state.  This  technique  was  introduced  by  the  Syage 
and  Wessel  in  1987  when  they  studied  the  dissociation  spectroscopy  of  napthalene 
ions  [163].  The  Ito  group  in  Japan  demonstrated  the  technique  subsequently  in  the 
spectroscopy  of  benzene,  p— difluorobenzene  and  1,3,5-trifluobenzene  [170].  Schlag’s 
group  demonstrated  ion  dissociation  spectroscopy  in  CH3I+,  CD3I+,  benzene,  fluo- 
robenzene,  [164,  165,  171,  172],  Since  then  a  number  of  different  groups  have  applied 
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this  technique  with  great  resolution  on  the  spectroscopy  [58,  61,  63,  64,  81,  166- 
168,  173-178], 

3. 7  Theory  of  two-photon  dissociation 

This  formalism  is  based  on  [179-182],  The  dynamics  of  two-photon  dissociation  can 
be  described  by  the  evolution  of  the  molecular  system  upon  interaction  with  laser 
radiation.  We  begin  with  the  time-dependent  Schrodinger  equation, 

ih—  =  TTT.  (67) 

dt  y  1 

Where  H  =  H0  +  Hi  and  H0  is  the  radiation  free  Hamiltonian  of  the  molecule, 

while  Hi  represents  the  interaction  with  the  radiation.  Following  the  treatment  of 
the  earlier  section,  now  we  have  two  lasers  coupling  three  states.  We  have  a  pulse 
laser  described  by  the  electric  field,  e\  =  e  ■  e(t)  cos(uot)  and  a  CW  laser  described  by 
62  =  eocos(a;if),  where  e  is  the  polarization  of  the  pulse  laser,  e(t)  is  the  slow  varying 
amplitude  of  the  electric  held  of  the  pulse  laser,  eo  is  the  amplitude  of  the  electric 
held  CW  laser.  The  interaction  Hamiltonian  follows  as, 

Hi  =  -fi0-e0(t)'-fii-e1(t).  (68) 

The  wavefunction  describing  the  process  becomes  a  linear  combination  of  the  radiation 
free  wavefunctions,  which  is  given  by, 

=  c0\^0)e-lEoh  +  ci\^i)e~iE^  +  J2  [  dE  cE,n\E,n-)e-iEt'\  (69) 

n 

where  the  coefficients,  c0,Ci,  and  cE,n  are  time  dependent  and  satisfy  the  condition 
that  (vh(f)|T(f))  =  1.  We  insert  the  wavefunction  into  the  time-dependent  Schrodinger 
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equation  67  to  get  a  set  of  first  order  differential  equations  for  the  coefficients  as, 

-^c0  =  in*0e~lAotCl(t),  (70) 

^-ci  =  in0e~lAotc0(t)  +i  I  dE^2nhEtne^AEtcE!n(t),  (71) 

n 

jtcE.n  =  i^lE,ne-iAEtCl(t), 

where, 

n  (^il/iol^o)e0(t)  n  m  ('0i|/^1|^7,  n-)ei(t) 

“0 \t)  —  - Z - ,  —  - Z - 

n  n 

Ei  —  E0  E  —  Ei 

A0  — - z - u)  o,  Ai  —  - - - u\.  (74) 

n  a 

We  can  solve  for  the  continuum  coefficient  by  imposing  the  boundary  condition  that 
prior  to  the  lasers  being  turned  on,  there  is  no  population  in  the  excited  state  and 
the  continuum.  We  obtain, 

CE,n(t)  =  i  [  dt'  ^l,E,n(t')elAEt'Cl(t')-  (75) 

Jo 

We  can  eliminate  the  continuum  coefficient  from  equation  71  by  substituting  the 
solution  from  equation  75  to  obtain, 

jcx  =  itt(t)eiAo  -  j  (IE  ]T  VthE,ne-lAEt  x  J*  dt'^E^(t')eiAEt' a(t').  (76) 

To  solve  the  above  equation,  we  make  the  assumption  that  all  the  matrix  elements  of 
the  continuum  states  vary  slowly  with  energy  i.e.  the  states  are  “flat”.  This  means 
that  energy  of  each  continuum  state  can  be  replaced  by  some  average  energy  such 
that  El  =  Ei  +  huji, 

(EL,n~\di\^i)\2  (77) 

n  n 

This  is  assumption  is  called  the  slowly  varying  continuum  approximation  (SVGA). 
At  this  point  it  becomes  intuitive  to  express  the  equation  in  terms  of  the  spectral 


(72) 


(73) 
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autocorrelation  function,  F\  (t  —  t '),  which  is  defined  as  the  Fourier  transform  of  the 
absorption  spectrum, 


Fx{t  -t')  =  J  dE  Ai(JS)e-iWB'l(t-t'),  (78) 

where  A\(E)  is  the  absorption  coefficient  from  the  ith  state  expressed  as, 

Al(E)  =  YJ\(EL^\^i\A)\2-  (79) 

n 

The  SVGA  localizes  the  autocorrelation  function  such  that  defining  it  in  terms  of 
equations  77  and  78,  we  find, 

Fi(t  -  if)  =  2'KhAl(EL)5{t  -  t').  (80) 


If  we  rewrite  the  differential  equation  for  c i  in  terms  of  the  autocorrelation  function, 
and  then  integrate  the  function  over  E  and  t',  we  obtain, 

-^ci  =  in,0(t)ezAotc0(t)  -  fti(t)ci(t),  (81) 


where  Oi  (t)  is  the  Rabi  frequency  given  by, 

n  (A;L,n_|^i|V'i)eiW 

“no  =  n2^ - 1 - 

n 


(82) 


The  probability  of  observing  the  dissociation  into  state  |n)  over  some  fixed  energy  E, 
becomes  the  square  of  the  photodissociation  coefficient, 

^  f  dt'(Oe^KAB*'ci(t') 

Jo 

In  our  REMPD  experiments,  we  observe  the  change  in  the  probability  of  dissociating 
the  molecule  as  a  function  of  the  detuning  from  the  vibrational  transition. 


Pn(E)  =  \cn>E(t  -)•  oo)|2  = 


Wi\Vi\E 


h 
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CHAPTER  IV 


EXPERIMENTAL  SETUP 


In  this  section  the  instruments  used  for  the  experiments  will  be  discussed  in  detail 
including  the  lasers  systems  used  and  the  control  of  the  laser  frequencies. 

4-1  Calcium  ion  diode  laser  system 

Diode  lasers  are  an  important  tool  in  atomic  physics.  Due  to  their  low  costs  compared 
to  dye  lasers  and  good  stability,  diode  lasers  are  lasers  of  choice  for  most  atomic 
physicists.  Briefly,  a  diode  laser  operates  by  emitting  photons  when  current  is  run 
through  the  active  region  of  between  a  n-type  and  p-type  cladding  layers.  This  creates 
electron-hole  pairs  that  recombine  to  emit  photons.  The  energy  of  the  emitted  photon 
is  determined  by  the  band-gap  of  the  semiconductor  [183].  Typical  semiconductors 
used  are  AlGaAs  for  the  750-890  nm  region,  InGaAsP  for  the  1300-1500  nm  region 
and  GaN  for  violet  light  with  a  range  365-405  nm. 

For  Doppler  cooling  and  sideband  cooling  40Ca+,  continuous  wave  (GW)  lasers 
are  used.  The  lasers  used  are  diode  lasers  (Toptica  DL100),  which  are  advantageous 
for  research  because  of  their  ease  of  use,  stability,  and  relatively  low  cost.  While  the 
design  of  the  lasers  is  the  same,  the  use  of  different  semi-conductors  gives  a  range  of 
available  wavelengths.  For  Toptica,  lasers  operating  wavelengths  with  a  range  370- 
1770  nm  are  available.  For  the  main  cooling  laser  for  40Ca+  at  397  nm  a  second 
harmonic  generated  frequency  multiplied  laser  is  used  with  the  fundamental  lasing  at 
794  nm.  A  typical  Toptica  laser  is  shown  in  the  picture  below  in  Fig.  18. 
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Figure  18:  A  picture  of  the  inside  of  a  Toptica  laser  (DL  100).  The  blue  arrow  shows 
the  direction  of  the  laser  beam  coming  out  of  the  laser.  The  wavelength  is  coarsely 
adjusted  with  the  screw  and  the  grating  position.  The  feedback  screws  are  used  to 
optimize  lasing. 


The  lasers  are  typically  coupled  to  the  traps  with  optical  fibers  from  another  opti¬ 
cal  table  since  the  lasers  are  shared  between  experiments.  Because  the  main  Doppler 
cooling  transition  requires  greater  precision  for  the  experiments,  an  additional  in¬ 
strument  is  used  to  control  the  frequency  of  the  397  nm  laser,  which  also  has  the 
ability  to  turn  on  and  off  the  laser.  For  our  experiments  we  used  an  acousto-optic- 
modulator  (AOM,  Brimrose  7998).  An  AOM  can  control  the  frequency,  intensity,  and 
the  direction  of  the  laser  beam.  The  device  works  by  Bragg  diffracting  incoming  light 
with  acoustic  wavepackets  that  are  propagating  through  the  crystal.  The  amplitude 
and  frequency  of  the  acoustic  wavefronts  can  be  changed  to  modulate  the  diffracted 
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light.  For  more  information  on  the  operation  of  AOMs,  the  reader  is  directed  to 
Ref.  [184], 

For  our  experiments  we  used  a  double  pass  configuration  to  set  up  the  AOM. 
Since  the  beams  are  fiber  coupled  to  the  traps,  any  change  in  the  angle  of  the  beam 
entrance  to  the  fiber  coupler  results  in  loss  of  intensity  on  the  trap  side.  However, 
scanning  the  frequency  of  the  AOM  results  in  changes  of  the  angle.  To  overcome 
this,  a  double  pass  configuration  is  utilized.  In  experimental  setups,  the  first  order 
positive  diffracted  light  which  is  retro-reflected  through  the  AOM  for  a  double  pass 
is  used.  This  removes  the  angle  shift  in  the  output  beam  such  that  the  input  beam  is 
co-aligned  with  the  output  beam.  To  separate  the  beams  we  use  polarization  optics 
as  shown  in  Fig.  19.  Light  is  first  passed  through  a  half-wave  plate  which  rotates 
the  beam  by  45  degrees.  It  then  passes  a  polarizing  beam  splitter  cube  before  a  lens 
focuses  it  into  the  AOM.  The  diffracted  beam  on  the  other  side  of  the  AOM  is  then 
re-collimated  with  a  second  lens.  The  wanted  diffraction  order  is  selected  with  an 
iris  while  the  rest  are  blocked.  That  light  passes  through  a  quarter-wave  plate  twice, 
which  makes  the  output  beam  orthogonally  polarized  to  the  incoming  beam. 
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Figure  19:  A  schematic  showing  the  optical  setup  for  the  double  pass  aucosto-optic 
modulator  (AOM).  The  setup  uses  polarization  to  separate  input  and  output  beams 
which  are  co-aligned.  A/2  is  a  half-wave  plate,  PBS  is  a  polarizing  beam  splitter 
cube, and  A/4  is  a  quarter- wave  plate.  The  iris  is  used  to  block  unwanted  beams  from 
being  retro-reflected 

4.1.1  Wavemeter 

The  frequency  of  the  CW  lasers  is  monitored  with  a  HighFinesse  WS/7  wavelength 
meter,  which  is  housed  in  a  thermally  isolated  casing.  The  instrument  has  an  absolute 
frequency  accuracy  of  60  MHz  and  a  frequency  range  of  350-1120  nm  and  acquisition 
speed  of  150  Hz.  However,  since  we  have  to  monitor  many  lasers  at  once  during  the 
experiments,  an  external  fast  operating  switch  is  coupled  via  a  multi-mode  fiber  to 
the  wavemeter.  This  reduces  the  the  absolute  accuracy  of  the  wavemeter  to  a  relative 
accuracy  of  200  MHz.  Once  a  month  or  so  the  wavemeter  would  drift  due  to  thermal 
fluctuations  of  the  room  and  the  wavemeter  was  re-calibrated  with  a  HeNe  laser  while 
an  ion  was  trapped.  Typical  measured  frequency  changes  after  recalibrations  were  on 
the  order  of  100s  of  MHz. 
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Jf.,2  Femtosecond  laser  system 

For  the  resonance  enhanced  mnltiphoton  dissociation  experiments,  we  required  a 
broadband  laser  that  could  address  all  the  rotational  levels.  We  use  a  mode-locked 
femtosecond  Titanium-Sapphire  (TkSapph)  laser.  For  our  experiments,  we  used  a 
Coherent  Mira  Optima  900-F  femtosecond  laser.  The  laser  is  pumped  by  a  Verdi- 
V10  with  15  W  of  power  and  outputs  about  2  W  of  mode-locked  power  that  can  be 
scanned  from  700  nm  to  980  nm.  Pulses  from  the  lasers  come  out  at  a  repetition  rate 
of  76  MHz  and  have  a  length  below  200  fs  with  peak  power  of  176  kW. 

To  check  that  the  pulse  is  not  distorted  when  it  passes  through  optics,  we  measured 
the  characteristics  of  the  laser  pulses  using  frequency-resolved  optical  grating  (FROG) 
method  provided  to  us  by  Rick  Trebino’s  group  [185,  186].  The  beam  was  measured 
at  different  locations  before  and  after  the  trap.  Briefly,  a  pulsed  laser  has  an  intensity 
and  phase  vs  time  (space),  which  can  be  described  with  an  electric  held  written  in 
the  time  domain  as, 

E(t)  =  Re  [  \/I{i)eiuJot-^] ,  (84) 


where  J(t)  is  the  intensity  and  0(f)  is  the  phase.  In  the  frequency  domain,  the  held 
can  be  written  as, 

E(u)  =  y/Siuje-^,  (85) 


where  S(u)  is  the  spectrum  and  (p(w)  is  the  spectral  phase.  The  FROG  apparatus 
measures  both  shaped  and  unshaped  pulses  by  measuring  the  intensity  and  phase  of 
a  femtosecond  pulses.  The  FROG  trace  can  be  characterized  as, 

r  oo  -i  2 


Ifrog{u,t) 


E(t)Eg(t  -  T)e~iujtdt 


(86) 


OO 


where  E(t)  is  the  electric  held  of  the  unknown  pulse,  u  is  the  frequency,  and  r  is 
the  delay.  Thus  a  FROG  measurement  shows  spectral  resolved  autocorrelations  of 
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intensity  versus  the  frequency  and  the  delay  to  characterize  the  pulse.  A  typical 
FROG  trace  is  shown  Fig.  20. 


Figure  20:  A  FROG  trace  of  the  laser 


The  data  from  the  FROG  measurements  is  shown  in  the  two  plots  below  in  Fig.  21. 
The  temporal  spectrum  in  Fig.  21a  shows  the  intensity  as  a  function  of  the  delay, 
which  gives  information  about  the  pulse  length.  According  to  the  plot  the  pulse 
length  is  163.2  fs,  which  is  in  agreement  with  the  manufacture  specifications.  The 
frequency  spectrum  in  Fig.  21b  shows  the  width  of  the  peak  of  the  laser  as  a  function 
of  wavelength.  At  this  central  frequency,  the  width  was  8.51  nm. 
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Figure  21:  Spectrum  from  the  FROG  profile,  (a)  shows  the  autocorrelation  between 
the  intensity  and  the  delay  while  (b)  shows  the  autocorrelation  between  the  intensity 
and  the  frequency. 
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To  measure  the  frequency  of  the  laser  during  the  experiment,  we  used  Ocean- 
Optics  HR2000+CG  spectrometer,  which  shows  only  the  intensity  as  a  function  of 
the  frequency. 

4-3  Vacuum  chambers 

In  order  to  trap  and  laser-cool  ions,  we  need  to  put  them  in  an  isolated  environment 
where  collisions  are  minimized.  For  this  reason,  ultra-high  vacuum(UHV)  chambers 
are  used  in  tandem  with  RF  Paul  traps.  Typically,  the  stainless  steel  parts  connected 
together  with  OF  sealing  are  used  to  build  the  chamber.  To  maintain  low  pressures, 
an  ion  pump  is  used  (Duniway  Stockroom/Gamma  Vacuum)  with  pump  speeds  of 
50  L/s  or  more.  An  additional  titanium  sublimation  pump  is  also  used  to  remove 
any  excess  particles  inside  the  chamber,  which  helps  maintain  low  pressures  inside 
the  chamber.  In  the  Brown  lab,  two  UHV  chambers  with  two  different  traps  were 
used  for  molecular  experiments.  For  historical  reasons,  we  called  them  the  molecular 
chamber  and  the  atomic  chamber.  A  drawing  of  the  chamber  housing  each  trap  is 
shown  in  Fig.  22. 
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Figure  22:  A  schematic  showing  an  assembled  chamber.  The  arrows  represents  the 
components  used  to  trap  ions  such  as  the  RF  and  DC  feedthroughs,  the  lasers,  the 
vacuum  system,  and  the  H2  leak  valve. 

4.3.1  Molecular  Chamber 

Laser-cooled  atomic  ions  as  probes  for  molecular  ions 

Kenneth  Brown,  C.  Ricardo  Viteri,  Craig  R.  Clark,  James  E.  Goeders,Ncamiso  B. 
Khanyile,  and  Grahame  D.  Vittorini,  AIP  Conf.  Proc.  1642,  392  (2015) 

We  setup  and  initially  trapped  in  the  Molecular  Trap  which  had  an  eleven  segment 
trap  designed  by  one  of  the  former  graduate  student,  James  Goeders  [187].  The  trap 
is  a  11-segmented  trap  which  makes  it  possible  to  trap  in  five  different  regions.  The 
electrodes  are  Nitronic  50  stainless  steel  with  curvatures  of  about  0.5  mm.  The  trap 
has  a  radius  (r0=1.0  mm)  where  ions  are  trapped.  Trapping  RF  voltage  is  applied  to 
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diagonally  opposed  electrodes  while  the  other  electrodes  are  kept  to  RF  ground.  To 
confine  the  ions  axially,  DC  voltage  is  applied  to  nine  3-mm-long  electrodes  diagonally 
opposite  from  each  other  while  there  are  additional  two  10-mm-long  endcaps.  To  trap 
ions,  the  trap  was  driven  at  0/2=10.41  MHz  with  400V  of  RF  voltage. 


6.35  mm 


Figure  23:  A  schematic  showing  the  11  segmented  Paul  trap,  we  call  the  Molecular 
Trap,  (a)  Shows  the  cross  section  of  the  electrodes  while  (b)  shows  a  3D  representation 
of  trap. 

4.3.2  Atomic  Chamber 

All  the  other  experiments  described  in  this  thesis  were  done  in  the  other  chamber, 
which  we  historically  called  the  Atomic  Chamber.  The  experiments  include  trap¬ 
ping  Yb+  ions  with  the  different  isotopes  (Chapter  2),  sympathetic  heating  spec¬ 
troscopy  [188],  sideband  cooling  of  40Ca+  and  40CaH+  [189],  and  measuring  the  vi¬ 
brational  overtones  of  10CaH+  (Chapter  5).  The  trap  was  designed  by  the  Urabe 
group  in  Japan  [190].  The  trap  is  a  5-segmented  trap  with  a  radius  (ro=0.6  mm).  Two 
diagonally  opposed  non-segmented  stainless  steel  1-nnn  thick,  3-nnn  long  electrodes 
mounted  on  insulation  spacers  are  used  to  apply  RF  voltage.  The  plates  have  a 
cross-section  curvature  of  0.5  mm.  The  DC  voltage  is  applied  to  two  of  3-mm-long 
segmented  electrodes  that  are  arranged  at  intervals  of  0.5  mm  to  confine  the  ions 
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radially.  To  compensate  for  stray  electric  fields,  additional  voltage  is  applied  to  non- 
segmented  electrodes  and  the  electrodes  at  the  center  of  the  trap.  To  trap  40Ca+,  we 
drove  the  trap  with  200V  of  RF  voltage  at  12/2=14.34  MHz. 


(a)  RF  electrodes 


(b)  Full  trap 


Figure  24:  A  schematic  of  the  5-segmented  Urabe  trap,  (a)  shows  a  cross-section  of 
the  RF  electrodes  while  (b)  shows  a  3D  view  of  the  trap. 


4-4  Radio  frequency  drive 


To  trap  ions,  about  200  V  of  RF  is  applied  to  the  electrodes.  To  generate  the  needed 
voltage,  a  functional  generator  is  used  to  produce  a  sinusoidal  wave  at  the  desired 
frequency.  The  produced  RF  is  the  amplified  with  a  commercial  amplifier  before  being 
introduced  into  the  trap.  To  introduce  the  voltage  in  the  trap,  a  helical  resonator  is 
used,  which  matches  the  impedance  between  the  trap  and  the  function  generator; 
enabling  high  voltages  to  be  introduced  in  the  trap.  In  addition,  the  resonator  also 
acts  as  a  filter  that  limits  any  noise  injected  into  the  trap  [191].  The  performance 
of  a  helical  resonator  is  determined  by  the  resonant  frequency,  lo0  =  and  the 
quality  factor,  Q  =  where  L  is  the  inductance  and  C  is  the  capacitance. 

When  designing  the  helical  resonator,  the  various  parameters  shown  in  Fig.  25  are 
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taken  into  account. 


Figure  25:  A  schematic  showing  the  design  of  a  helical  resonator.  Here  the  case 
diameter  is  shown  as  D,  the  case  height  as  B,  the  coil  diameter  d,  coil  height  b,  the 
winding  pitch  r,  the  coil  wire  diameter  do. 

To  design  the  helical  resonator  the  following  equations  are  used  with  the  param¬ 
eters  above.  In  terms  of  the  parameters,  the  quality  factor  can  be  express  as, 

Q  =  35 .9dy/f,  (87) 

where  /  is  the  resonant  frequency  of  the  resonator.  Other  characteristics  can  be 
calculated  with  the  following  equations, 
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where  N  is  the  number  of  turns  and  Z$  is  the  characteristic  impedance  of  an  air-filled 
helical  transmission.  The  characteristic  capacitance  resonator  is  given  by 

Cs  ~  bKCs(d,  D),  (91) 


where 


0  75 

Kn,  =  39.37,  —  x  10~12F/m. 


l°9h 


(92) 


Generally,  these  parameters  act  as  a  guide  for  the  design.  The  resonator  is  nor¬ 
mally  optimized  empirically  to  suit  the  needs  of  the  experiment.  For  the  resonator 
used  in  the  atomic  trap,  the  resonant  frequency  was  14.14  MHz  while  the  quality 
factor  Q  was  120.  The  Q  factor  is  relatively  high  for  a  resonator  at  this  low  frequency. 


4-5  Ion  detection 

Trapped  ions  are  detected  via  fluorescence  from  the  spontaneously  emitted  photons 
on  the  cooling  transition.  Since  the  scattered  photons  from  the  ions  are  scattered  in 
random  directions,  efficient  imaging  optics  with  high  numeral  apertures  (N.A.)  are 
needed  to  collect  the  photons  and  image  them  [116].  In  addition,  the  imaging  optics 
need  to  work  close  to  the  diffraction  limit.  The  collected  photons  are  then  imaged 
with  an  EMCCD  camera  and  counted  with  a  photon  multiplier  tube  (PMT).  For  our 
experiments  the  lens  system  was  designed  and  optimized  with  Zeemax  commercial 
software.  A  schematic  of  the  lens  system  used  to  collect  photons  is  shown  in  Fig.  26. 
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Figure  26:  A  design  of  the  lens  system  to  image  fluorescence  from  a  4.5  inch  octagon. 
The  lenses  used  are  NBK-7  coated  to  minimize  back  reflection  from  the  lens  windows 

The  lens  system  had  an  N. A. =0.43,  a  magnification  of  lOx  and  had  a  total 
length  of  660  mm.  The  lenses  used  were  bought  from  Thorlabs  inc,  and  they  were 
LAl401(d=2  inches,  /= 60.0  mm),  LC1611  (d= 2  inches,  /=-150.0  mm),  LA1384  (d= 2 
inches,  /= 125.0  mm),  LA1399  (d= 2  inches,  /= 175.0  mm),  where  d  is  the  diameter 
of  the  lens  and  /  is  its  focal  length.  The  lenses  are  listed  in  their  order  of  appearance 
from  the  trap. 

4.6  Frequency  stabilization  of  laser  systems 

For  laser  cooling  experiments,  the  frequency  of  the  laser  must  not  drift  by  more  than 
the  linewidth  of  the  transition  (^20  MHz).  However,  the  laser  frequency  is  affected 
by  environmental  factors  such  as  the  fluctuations  in  temperature  and  pressure.  To 
counter  this,  an  external  cavity  is  used  to  lock  the  frequency  of  the  laser. 

Lasers  are  often  locked  to  one  of  the  modes  of  a  Fabry-Perot  interferometer 
(FPI)  [192],  However,  because  an  FPI  cavity  needs  to  be  actively  temperature  con¬ 
trolled  and  under  vacuum  to  minimize  pressure  fluctuations,  it  makes  the  setup  expen¬ 
sive.  A  cheaper  alternative  is  to  use  a  cavity  in  ambient  pressure,  which  is  referenced 
to  a  stable  laser  [193,  194].  In  our  lab,  we  use  a  stabilized  HeNe  laser  with  a  drift  of 
about  2  MHz/8  hrs,  which  is  ideal  for  laser  cooling  experiments. 

The  transfer  cavity  works  by  aligning  a  stable  laser,  in  our  case  a  HeNe,  and 
unstable  laser(s)  to  a  cavity.  Light  transmitted  from  the  cavity  is  monitored  while  the 
cavity  length  is  scanned  with  a  piezo  actuator,  which  produces  peaks  as  a  function 
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of  the  distance  traveled  by  the  piezo  [195].  The  cavity  is  locked  to  the  position  of  the 
HeNe  peak  while  the  lasers  arc  locked  to  their  position  through  feedback  to  the  laser. 

Light  enters  the  cavity  by  leaking  through  a  backside-polished  mirror.  The  light  is 
reflected  between  the  two  cavity  mirrors.  As  the  piezo  is  scanned,  light  in  the  cavity 
is  amplified  when  the  mirror  spacing  is  in  resonance  with  an  integer  number  of  the 
half- wavelengths  like, 

2L/X  =  n,  (93) 

where  L  is  the  length  of  the  cavity,  A  is  the  wavelength,  and  n  is  an  integer.  This 
causes  constructive  interference  when  Sn  =  ±1.  The  intensity  of  the  light  inside 
the  cavity  increases  and  leaks  through  the  back  mirror  where  it  is  detected  with  a 
photodetector. 

4.6.1  Realization  of  a  transfer  cavity 

We  locked  three  lasers, 397  nm,  866  nm,  and  854  nm  with  a  transfer  cavity  to  the 
HeNe.  Since  the  397  nm  laser  is  generated  by  frequency,  we  used  the  fundamental  to 
lock  it  because  the  IR  wavelengths  have  better  finesse  with  the  mirrors  we  used.  The 
setup  and  the  cavity  are  shown  below  in  Fig.  27. 
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Figure  27:  A  schematic  showing  the  optical  layout  for  the  transfer  cavity.  Lasers  are 
coupled  into  the  cavity  using  a  combination  of  polarization  and  band-pass  Liters.  An 
achromatic  lens  in  front  of  the  cavity  mode-matches  the  lasers  to  the  cavity.  The  light 
transmitted  from  the  cavity  is  sent  to  detectors  for  signal  processing  and  feedback  to 
the  lasers. 


An  important  feature  of  coupling  light  into  the  cavity  is  the  use  of  polarization 
and  band  pass  Liters  to  combine  and  separate  the  diLerent  wavelengths.  Since  the 
794  nm  and  the  633  nm  of  the  laser  are  relatively  close  in  wavelength,  we  use  a 
polarization  beam  splitter  (PBS)  cube  to  combine  them  and  we  do  the  same  for  the 
854  nm  and  the  866  nm.  To  combine  the  854/866  nm  and  the  794/633  nm  lasers,  we 
use  a  dichroic  mirror  (Thorlabs  DMSP-805).  The  lasers  are  then  mode-matched  to 
the  cavity  using  an  achromatic  lens  where  the  lasers  are  focused  to  the  back  mirror. 
The  cavity  is  a  half-concentric  etalon  with  the  Lrst  mirror  being  a  backside  polished 
concave  mirror  (Thorlabs  CM254-075-E02P)  and  the  back  mirror  is  a  Lat  broadband 
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(Semrock  MaxMirror)  mirror  that  has  over  99%  reflectivity  over  a  wide  range  as 
shown  in  Fig.  28.  The  lasers  are  separated  in  the  same  manner  they  were  combined 
to  individual  detectors.  To  ensure  no  leakage  of  other  light  to  the  detectors,  narrow 
band  filters  are  placed  in  front  of  each  detector. 


Figure  28:  A  schematic  showing  a  cross  section  of  the  cavity.  The  cavity  is  made  of 
a  concave  mirror  and  a  flat  mirror  inside  a  1-inch  tube.  The  mirrors  are  separated 
by  a  low-thermal-expansion  quartz  spacer.  The  flat  mirror  is  scanned  with  a  piezo  to 
cover  the  free  spectral  range  of  the  lasers. 

The  cavity  is  housed  in  a  one  inch  tube  and  the  two  mirrors  are  separated  with 
a  low-expansion  quartz  spacer  ( 8L/L  «  0.54  x  10-6  K-1),  three  piezos  are  mounted 
on  the  back,  flat  mirror  and  a  sinusoidal  voltage  of  50  V  is  applied  at  30  Hz.  This 
gives  the  cavity  a  scanning  range  of  about  2.1  pm,  which  is  enough  to  cover  the  free 
spectral  range  of  all  the  lasers.  A  scan  showing  the  transmission  of  the  397  nm  laser 
as  the  piezo  is  scanned  is  shown  in  Fig.  29. 
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Figure  29:  Transmission  of  the  397  nm  laser  through  the  cavity.  The  peaks  occur 
when  the  piezo  is  in  a  position  where  it  is  on  resonance  with  half-integer  wavelengths 
of  the  laser.  The  distance  between  two  peaks,  AA,  is  called  the  free  spectral  range. 
The  full  width  half  maximum,  5 A,  is  the  linewidth  of  the  cavity.  Together  these  two 
parameters  describe  the  finesse  of  the  cavity  by  /  =  where  R  is  the 

reflectivity  of  the  mirrors.  For  our  transfer  cavity,  the  finesse  is  «120. 


4.6.2  Locking  performance 

To  lock  the  lasers  we  used  home-made  photodiode  detectors,  along  with  National 
Instruments  cards  to  collect  the  data,  process  the  error  signal  and  send  feedback  to 
the  laser.  To  evaluate  the  performance  of  the  cavity,  we  looked  at  the  frequency  the 
397  nm  laser  as  a  function  of  time  is  shown  below  in  Fig.  30. 
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Figure  30:  A  plot  showing  the  frequency  of  a  locked  397  nm  laser  over  time.  The 
spikes  in  the  frequency  are  from  sudden  pressure  changes  in  the  room  such  as  someone 
opening  and  closing  the  door.  The  black  trace  is  a  running  average  of  the  individual 
blue  data  points.  The  measurement  in  the  case  is  also  limited  by  the  accuracy  of  the 
multi-channel  switch  of  the  wavemeter,  which  has  a  relative  accuracy  of  ±200  MHz. 


To  test  the  stability  of  the  locked  laser  we  plot  the  Allan  variance  of  the  frequency 
of  the  laser  over  time.  The  Allan  variance  is  an  averaging  statistic  used  to  measure  the 
stability  of  an  oscillator  [196,  197].  The  Allan  variance  is  described  by  the  equation, 


ay{T)  = 


M- 1 


XA+i  -  2*)2- 


(94) 


\|  2<M  -  !)  ± 

where  y%  describes  frequency  offset1  measurements  of  the  yi,  y2,  ■■■Vn  values,  M  is  the 


lrThis  is  the  deviation  of  the  measured  frequency  from  the  ideal  frequency. 
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number  of  values  in  the  series,  and  r  is  the  segment  between  the  data  points  measured 
in  seconds.  A  plot  of  the  Allan  variance  shows  the  improvement  of  the  measurement 
as  the  averaging  period  (r)  is  increased.  To  test  the  stability  transfer  lock  we  input  the 
measurements  of  the  397  nm  laser  frequency  from  the  WS7  in  to  the  Allan  equation 
and  we  get  the  plot  shown  in  Fig.  31. 
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Figure  31:  An  Allan  variance  of  the  locked  397  nm  laser,  which  was  locked  to  the 
transfer  cavity.  The  plot  shows  that  the  deviation  in  the  laser  reaches  about  2  MHz 
over  the  course  of  a  day.  The  slope  of  the  Allan  variance  identifies  the  type  of  noise  in 
the  system.  r_1  is  flicker  noise  and  white  noise,  r-1/2  is  the  white  frequency  noise  and 
r°  is  the  flicker  frequency.  The  Allan  variance  shows  the  removal  of  flicker  noise  and 
white  noise  as  we  increase  the  averaging.  Longer  averaging  though  can  increase  the 
noise  in  the  form  of  random  walk  noise.  In  this  graph,  if  we  continued  to  average,  we 
would  expect  the  slope  to  be  positive.  Normally,  we  would  expect  that  increasing  the 
averaging  time  further  than  tail0  would  result  in  poor  performance,  hence  the  graph 
would  start  going  up.  In  our  example,  it  looked  like  we  had  stopped  the  experiment 
prematurely. 


We  can  see  from  these  plots  that  the  transfer  cavity  effectively  stabilizes  the  laser 
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over  a  long  time.  The  spikes  in  the  frequency  data  represent  short  term  perturbations 
such  as  an  opening  and  closing  of  the  door,  which  the  cavity  does  not  respond  fast 
enough  to.  However,  during  the  course  of  a  day  we  see  stabilities  of  ±4  MHz,  which 
is  sufficient  for  laser  cooling  experiments. 
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CHAPTER  V 


OBSERVATION  OF  VIBRATIONAL  OVERTONES  BY 
SINGLE  MOLECULE  RESONANT 
PHOTODISSOCIATION 


This  chapter  is  currently  in  review. 

Ncamiso  Khanyile,  Gang  Shu  and  Kenneth  Brown,  in  review  (2015) 

5. 1  Introduction 

Precision  spectroscopy  of  molecules  and  molecular  ions  can  yield  insight  into  the  fun¬ 
damental  physical  constants  and  astrochemical  processes  [198,  199].  Coulombic  crys¬ 
tals  of  molecular  ions  and  laser-cooled  atomic  ions  are  ideal  systems  for  these  measure¬ 
ments  [98].  Cold  molecular  ions  have  been  used  to  study  cold  chemistry  [57,  103,  200], 
to  control  rovibrational  states  [66,  67,  70,  73,  201,  202],  and  to  precisely  measure 
molecular  transitions  [63,  203].  However,  expanding  these  techniques  to  a  wider  array 
of  molecular  ion  species  remains  a  challenge  clue  to  the  lack  of  experimental  data 
on  molecular  ions  transitions.  This  requires  new  methods  for  obtaining  spectral  in¬ 
formation.  The  spatial  localization  of  molecular  ions  in  a  Coulomb  crystal  results 
in  the  required  ion  density  for  spectroscopy  with  low  ion  numbers  relative  to  tradi¬ 
tional  techniques.  This  enhanced  sensitivity  has  recently  been  used  to  directly  observe 
dipole-forbidden  transitions,  whose  transition  frequencies  were  previously  inferred 
from  measurements  of  allowed  transitions  [203]. 

Coulombic  crystals  composed  of  laser-cooled  atomic  ions  and  molecular  ions  pro¬ 
vide  a  pristine  environment  for  studying  the  properties  of  molecules.  The  laser-cooled 
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atomic  ions  serve  as  both  a  coolant  that  reduces  the  temperature  and  a  sensitive  de¬ 
tector  that  allows  for  single  molecule  measurements.  The  system  has  been  used  to 
study  fundamental  reaction  dynamics  ranging  from  kinetic  isotope  effects  [103]  to  ob¬ 
serving  the  relative  reaction  rates  of  molecules  whose  shape  differs  by  the  orientation 
of  a  single  bond  [57].  It  is  also  a  natural  system  for  precision  measurements  of  molec¬ 
ular  ion  transitions.  Molecular  transitions  allow  for  the  precise  study  of  fundamental 
constants  using  physics  that  is  inaccessible  in  atomic  systems  [204,  205].  Spectroscopy 
with  Coulomb  crystals  has  been  used  to  directly  measure  dipole-forbidden  transitions 
in  Nj  [203]  and  accurately  determine  the  rovibrational  spectrum  of  HD+  [63].  Future 
experiments  to  test  the  stability  of  fundamental  constants  in  time  have  been  pro¬ 
posed  [28]  using  the  two-ion-species  techniques  already  exploited  in  the  most  precise 
atomic  ion  clocks  [31,  121], 

Systems  built  for  high-precision  measurement  are  often  incompatible  with  the 
survey  spectroscopy  required  to  find  unknown  transitions.  However,  this  is  not  the  case 
for  Coulomb  crystals,  where  the  long  ion  storage  time  provides  multiple  opportunities 
to  probe  the  molecule  and  the  fluorescence  of  the  laser-cooled  atomic  ion  serves  as 
a  fast,  low-noise  detector.  40CaH+  is  a  candidate  molecule  for  testing  the  possible 
variation  in  the  proton-to-electron  mass  ratio  [141] .  It  is  also  expected  to  be  relatively 
abundant  in  space  due  to  observation  of  CaH,  but  has  not  yet  been  directly  observed 
[17].  In  both  cases,  laboratory  measurements  of  bound  transitions  in  40CaH+  are 
required  for  scientific  progress 

Here  we  measure  two  previously  unobserved  vibrational  overtones  of  40CaH+  with¬ 
out  any  prior  experimental  data  to  guide  our  search.  We  use  small  crystals  of  40Ca+ 
and  40CaH+  to  observe  vibrational  transitions  via  resonance  enhanced  multiphoton 
dissociation  detected  by  40Ca+  fluorescence.  Based  on  theoretical  calculations  [25],  we 
assign  the  observed  peaks  to  the  transition  from  the  ground  vibrational  state,  v  =  0, 
to  v  —  9  and  v  =  10.  Our  method  allows  us  to  track  single  molecular  events,  and  it 
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can  be  extended  to  work  with  any  molecule  by  using  normal  mode  frequency  shifts  to 
detect  the  dissociation.  This  survey  spectroscopy  serves  as  a  bridge  to  the  precision 
spectroscopy  required  for  molecular  ion  control.  Our  results  are  an  important  step 
towards  using  molecular  vibrations  to  test  possible  time  variation  in  the  proton-to- 
electron  mass  ratio  [28,  29,  141]  and  to  control  cold  molecular  reactions  one  molecule 
at  a  time. 

5.2  Methods 

We  use  laser-cooled  40Ca+  ions  to  sympathetically  cool  40CaH+  molecules.  The  40Ca+ 
ions  are  trapped  in  a  RF  Paul  trap  (r0  =  0.6  mm)  driven  at  O  =  27r  x  14  MHz  to 
confine  the  ions  radially,  while  static  DC  voltage  applied  at  the  endcaps  confines  the 
ions  axially.  The  trap  is  kept  at  a  base  pressure  of  about  4  x  10~9  Pa.  The  ions 
are  detected  by  laser  induced  fluorescence  at  397  nm  onto  a  photon  multiplier  tube 
(PMT)  and  electron  multiplying  charge  coupled  device  (EMCCD)  camera.  A  narrow 
bandpass  filter  is  used  to  ensure  that  only  397  nm  light  is  detected. 

A  40CaH+  molecule  is  produced  by  leaking  about  5  x  10~7  Pa  of  molecular  H2  into 
the  chamber  via  a  leak  valve.  A  chain  consisting  of  three  40Ca+  ions  is  trapped  before 
leaking  in  the  H2.  The  40CaH+  is  produced  via  inelastic  collisions  in  the  gas  phase 
between  40Ca+(4Pi/2)  and  the  H2  as  40Ca+  +  H2  — »  40CaH+  +  H.  The  occurrence 
of  a  reaction  is  determined  when  one  of  the  ions  goes  dark  and  there  is  a  drop  in 
fluorescence  counts.  Once  a  reaction  occurs,  the  leak  valve  is  closed  and  the  experi¬ 
ment  is  delayed  until  the  base  pressure  is  reached.  After  this  pump  down  time,  the 
internal  degrees  of  freedom  of  the  molecule  are  expected  to  be  at  room  temperature 
due  to  weak  coupling  with  the  cold  translational  degrees  of  freedom.  The  identity  of 
the  molecule  can  be  determined  by  resolved  sideband  spectroscopy  [206]  and  under 
these  experimental  conditions  we  have  only  observed  the  formation  of  40CaH+. 

To  measure  the  vibrational  overtones  of  40CaH+,  we  use  (1  +  F)  REMPD.  A 
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broadband  mode-locked  Ti:Sapph  (IR)  laser  excites  a  specific  vibrational  transition 
of  40CaH+,  a  second  photon,  a  fixed  ultraviolet  (UV)  cw  laser  (380  nm)  dissociates 
the  ion  from  the  upper  vibrational  state  by  coupling  it  to  a  repulsive  electronic  state. 
Upon  dissociation,  the  previously  dark  40CaH+  ion,  will  be  broken  into  40Ca+  +  H, 
the  40Ca+  will  fluoresce  again  and  there  will  be  an  increase  in  fluorescence  counts 
as  shown  below  in  Fig.  33.  Other  possible  dissociation  channels,  including  two  UV 
photon  absorption  and  UV  induced  electron  bombardment,  are  measured  to  be  low 
compared  to  the  REMPD  dissociation  rates  by  performing  control  experiments  with 
the  IR  blocked.  A  pinhole  before  the  PMT  allows  partial  light  collection  from  all  three 
ions  and  reduces  background  due  to  scattered  light.  Misalignment  from  the  crystal 
center  results  in  three  distinct  collection  efficiencies  for  each  ion  position.  This  allows 
us  to  detect  the  position  of  the  dark  ion  from  the  fluorescence  as  shown  below  in 
Fig.  33. 
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Figure  32:  Energy  level  diagram  of  40CaH+.  Simplified  40CaH+  energy  level  diagram 
showing  the  overtones  excited  by  a  pulsed,  tunable  infrared  laser  (800-900  nm).  A  sec¬ 
ond  ultraviolet  laser  (377  nm)  excites  the  overtones  to  the  unbound  state  to  dissociate 
the  molecule. 

The  calculated  values  for  40CaH+  properties  are  from  Abe  et  al.  [25]  and  based 
on  the  method  NRel/cc-pCV5Z/CASPT2. 

Resonance  Enhanced  Multi-Photon  Dissociation  (REMPD)  is  a  common  tool  in 
physical  chemistry  where  conditional  on  a  resonant  transition  a  hnal  photon  dis¬ 
sociates  the  molecule.  It  can  be  applied  to  study  the  vibrations  of  molecules  from 
diatomics  with  a  single  electron  [64]  to  peptides  [207].  Ion  fragments  are  typically 
detected  by  an  electron  multiplier  after  mass  selection  [207,  208].  Embedding  the 
molecular  ions  in  laser-cooled  atomic  ions  opens  up  other  approaches  of  detection 
based  on  changes  in  the  Coulomb  crystal  shape  [209]  and  modulation  of  the  atomic 
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ion  fluorescence  by  driving  the  resonant  motion  of  the  molecular  ions  [64],  Resolved 
sideband  spectroscopy  on  a  narrow  transition  of  the  atomic  ions  can  also  be  used  to 
identify  the  fragments  [206].  For  the  case  of  40CaH+,  the  measurement  is  simplified 
by  the  trapped  fragment  40Ca+  fluorescing  brightly.  This  allows  us  to  use  the  ion  re¬ 
cycling  approach  that  was  previously  used  to  measure  the  reaction  kinetics  of  40Ca+ 
+  HD  by  direct  dissociation  [210]. 

Although  the  molecular  ion  is  at  a  translational  temperature  of  a  few  millikelvin, 
the  internal  degrees  of  freedom  are  in  equilibrium  with  the  room  temperature  vac¬ 
uum  chamber  via  black  body  radiation.  The  calculated  vibrational  frequency  of  the 
molecule  is  1478.4  cm-1  and  we  expect  that  the  molecule  will  be  in  the  ground  vi¬ 
brational  state  X4S+  greater  than  99.9%  of  the  time  [25].  On  the  other  hand,  the 
calculated  ground  state  rotational  constant,  4.711  cm-1,  is  small  relative  to  room 
temperature.  The  rotational  states  will  be  populated  with  an  expected  value  of  J  = 
5.36  and  the  lowest  ten  J  states  are  expected  to  have  more  than  94%  of  the  popula¬ 
tion.  Our  experiment  uses  a  single  molecule  at  a  time,  but  the  blackbody  radiation  will 
randomize  the  J  state  on  the  order  of  minutes.  To  cover  this  rotational  broadening  we 
use  a  150  fs  mode-locked  laser  to  drive  the  overtones.  Alternative  approaches  include 
reducing  the  rotational  state  space  by  sympathetic  cooling  with  a  buffer  gas  [70]  or 
a  cloud  of  laser-cooled  neutral  atoms  [73]  or  prepare  the  molecular  ions  at  a  specific 
state  by  threshold  ionization  [209].  Rotational  cooling  by  optical  pumping  [66,  67,  211] 
is  currently  impossible  for  40CaH+  due  to  the  lack  of  data  on  bound  transitions. 

The  experiment  takes  place  in  a  chamber  described  previously  [188].  We  start 
with  loading  three  40Ca+  ions  by  isotope  selective  photoionization  and  then  laser 
cooling  the  ions  to  form  a  three  ion  chain.  Hydrogen  gas  is  leaked  into  the  chamber  at 
5  x  1CT'  Pa  until  a  single  ion  reacts  as  detected  by  a  loss  of  fluorescence.  Occasional 
collisions  with  background  gas  will  rearrange  the  ion  chain  as  observed  by  the  relative 
position  of  the  dark  ion.  A  pinhole  in  front  of  the  photomultiplier  tube  allows  us  to 
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detect  these  collisions  by  modulation  in  the  photon  signal  as  shown  in  Fig.  33. 
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Figure  33:  Dissociation  measurement.  After  reaction  with  hydrogen,  a  mixed  Coulomb 
crystal  of  two  40Ca+  ions  and  one  40CaH+  ion.  The  dissociating  lasers  are  applied  and 
the  time  to  dissociation,  Td,  is  measured  by  observing  the  change  in  fluorescence.  The 
red  trace  shows  an  example  dissociation  event  when  the  infrared  laser  is  resonant  with 
an  overtone.  The  blue  trace  is  an  event  when  the  infrared  laser  is  blocked.  A  pinhole 
in  the  optical  system  allows  us  to  correlate  fluorescence  with  not  only  the  number  of 
40Ca+  ions  but  also  the  relative  position  of  the  molecular  ion. 

The  two  photon  photodissociation  occurs  by  first  exciting  a  vibrational  overtone  of 
X4S+  and  then  photodissociating  to  the  state  C4S  as  shown  in  Fig.  32.  We  drive  the 
overtone  using  a  femtosecond  mode-locked  spectroscopy  laser  with  5  nm  bandwidth, 
allowing  us  to  ignore  the  line  broadening  due  to  rotational  transitions.  The  second 
photon  comes  from  a  continuous- wave  laser  at  377  nm.  For  each  central  wavelength 
of  the  spectroscopy  laser,  we  expose  the  molecular  ion  to  the  dissociation  lasers  and 
record  the  time  for  the  molecule  to  dissociate,  Td  as  shown  in  Fig.  33.  The  transition 


80 


strength  is  proportional  to  1/  <  >. 


5.3  Results  and  discussion 


Figure  34:  40CaH+  vibrational  overtone  spectra.  The  measured  are  averaged  over 
eight  experiments  and  the  inverse  is  plotted  as  a  function  of  the  IR  wavelength.  Error 
bars  are  propagated  from  the  standard  deviation  of  (r^).  The  data  reveal  two  peaks 
which  are  fit  assuming  a  Gaussian  line  shape.  Gray  bars  are  centered  at  the  calculated 
theoretical  values [25]  for  the  z/  =  10  v  =  0  and  v’  —  9  ■<—  v  —  0  overtones. 

The  spectrum  clearly  shows  two  peaks  which  we  identify  as  the  v'  =  10  v  =  0 
and  v'  —  9  •<—  v  —  0  overtones  of  10CaH+  based  on  theoretical  calculations  [25]  as 
shown  in  Fig.  34.  We  measure  the  v'  —  10  ■<—  v  —  0  transition  to  be  at  812(3)  nm 
compared  to  the  theoretical  value  of  813.3  nm  and  v'  =  9  v  =  0  transition  to  be 
at  890(3)  nm  compared  to  the  theoretical  value  of  883.3  nm.  This  disagreement  is  in 
line  with  observed  differences  between  calculated  and  measured  vibrational  transition 
frequencies  in  other  metal  hydrides  [25]. 
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Our  experimental  setup  was  intended  for  high-precision  quantum  logic  spectroscopy 
[212]  experiments  on  molecular  ions.  We  have  shown  that  the  same  setup  can  be  used 
for  the  preliminary  large  range  spectroscopy  necessary  to  observe  even  weak  lines  de¬ 
spite  trapping  only  a  few  ions  at  a  time.  The  next  step  for  precision  spectroscopy  of 
40CaH+  is  to  reduce  the  rotational  temperature  by  sympathetic  cooling  with  neutral 
atoms  and  then  rotationally  resolve  these  transitions  and  the  fundamental  transi¬ 
tion  [70,  73].  Then  quantum  logic  spectroscopy  can  be  performed  on  ground  state 
cooled  40Ca+-40CaH+  crystals  [189]  in  order  to  reach  the  precision  necessary  for  ob¬ 
serving  relative  changes  in  fundamental  constants  [31]. 

Our  current  experiment  takes  advantage  of  the  fluorescence  of  the  dissociated 
product  and  can  be  easily  applied  to  other  fluorescent  product  ions.  The  method 
can  be  modified  for  dark  product  ions  by  periodically  blocking  the  dissociation  laser 
and  measuring  the  mass  of  the  product  ion.  This  can  be  achieved  by  monitoring  the 
sidebands  of  the  atomic  ions  and  calculating  the  mass  of  the  unknown  ion  from  the 
observed  normal  modes  [206].  We  expect  that  few  ion  Coulomb  crystals  will  provide 
a  method  for  discovering  the  spectra  of  astrochcmically  relevant  ions  that  are  difficult 
to  make  in  the  high  numbers  required  for  traditional  spectroscopic  methods. 

The  two  photon  photodissociation  occurs  by  first  exciting  a  vibrational  overtone  of 
X:S+  and  then  photodissociating  to  the  state  C4S  as  in  Fig.  32.  We  drive  the  overtone 
using  a  femtosecond  mode-locked  [pulsed]  spectroscopy  laser  with  5  nm  bandwidth, 
allowing  us  to  ignore  the  line  broadening  due  to  rotational  transitions.  The  second 
photon  comes  from  a  continuous-wave  laser  at  380  nm  For  each  central  wavelength  of 
the  spectroscopy  laser,  we  expose  the  molecular  and  record  the  time  for  the  molecule 
to  dissociate,  r^.  The  transition  strength  is  proportional  to  1/  <  >  Fig.  33. 
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5.4  Experimental  tests  and  controls 

There  are  different  mechanisms  that  can  lead  to  the  dissociation  of  the  molecule 
inside  the  trap.  To  ensure  that  the  signal  we  obtained  from  the  measurement  was 
due  to  resonant  dissociation  only,  we  conducted  some  experimental  controls.  The 
molecule  can  be  dissociated  by  background  collisions,  successive  excitation  caused  by 
absorbing  many  IR  photons,  electron  impact  dissociation  from  electrons  scattered  of 
the  electrodes  of  the  trap.  These  can  compete  with  REMPD  and  were  investigated. 

5.4.1  Infrared  multiple  photon  dissociation 

Prior  to  conducting  the  experiment,  we  worried  about  the  competition  between 
REMPD  and  infrared  multiple  photon  dissociation  (IRMPD).  This  is  a  process  where 
absorption  of  multiple  infrared  photons  can  induce  dissociation  of  strongly  bound 
molecules  [81,  156,  166,  176,  213].  IRMPD  requires  the  absorption  of  a  few  hundred 
photons;  hence,  very  strong  gas  discharge  lasers  are  used  such  as  CCD  lasers.  The 
MIRA  laser  has  very  high  peak  power  (176  Kilowatts)  and  delivers  about  1.13xlOn 
photons  per  pulse,  which  could  supply  the  necessary  number  of  photons  needed  for 
dissociation.  The  energy  needed  to  couple  the  ground  state  of  40CaH+  (X  1S+)  with 
the  repulsive  excited  electronic  state  (C  1II)  at  272  nm  could  be  achieved  by  three 
photon  excitation  close  to  the  v'  —  10  <—  v  —  0  transition  at  816  nm  as  shown  below 
in  Fig.  35a.  We  investigated  the  possibility  of  such  an  event  occurring  by  irradiating 
the  crystals  with  the  IR  laser  only  around  where  we  expect  the  peak.  The  data  was 
found  to  be  consistent  with  no  peak  as  shown  in  Fig.  35b. 
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(a)  Principle  of  IRMPD 
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Figure  35:  A  plot  of  IRMPD  data  and  principle,  (a)  shows  the  principle  of  infrared 
multiphoton  dissociation  while  (b)  shows  the  data  from  IRMPD  experiment.  The  data 
indicates  the  three-photon  process  is  indistinguishable  from  the  background  rate. 


Each  experiment  data  point  was  repeated  six  times  and  averaged.  The  resulting 
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spectra  show  an  essentially  flat  line  indicating  no  significant  triple  photon  dissociation 
of  the  the  molecular  ion.  In  addition,  the  dissociation  time  was  very  slow,  £s5xl02 
s.  Thus,  we  can  conclude  that  the  peak  was  not  caused  by  IR.MPD  but  rather  by 
REMPD.  This  could  be  because  the  excited  states  have  very  low  absorption  cross 
sections  and  there  might  not  be  virtual  states  present  for  multiple  photon  dissociation. 

5.4.2  Dissociative  recombination 

Initial  investigations  into  the  dissociation  spectra  were  conducted  with  two  pulsed 
lasers.  One  was  the  MIRA  while  the  other  one  was  the  Paladin  tripled  Nd:YAG 
laser  at  355  nm.  The  laser  combination  was  enough  to  couple  the  bound  state  to 
the  dissociative  state.  Initially,  only  the  1/  —  10  <—  v  —  0  vibrational  overtone  was 
investigated.  The  data  from  those  experiments  is  shown  below  in  Fig.  36.  The  data 
are  shown  in  the  red  circles  with  the  corresponding  error  bars.  The  blue  curve  is  a 
simulated  spectrum  given  the  width  of  the  laser.  The  data  and  the  simulated  spectra 
agree  very  well.  The  point  around  800  nm  differs  a  little  bit  from  the  spectra,  which  is 
similar  to  the  point  around  870  nm.  This  could  be  due  to  the  fact  that  the  simulation 
uses  an  average  width  for  the  laser  while  the  width  changes  as  the  wavelength  changes. 
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Figure  36:  Spectrum  of  the  v'  —  10  •<—  v  —  0  vibrational  overtone  obtained  by  using 
two  pulsed  lasers.  One  being  the  MIRA  and  for  dissociation,  a  tripled  Nd:YAG  at 
355  nm. 

We  switched  the  pulsed  laser  to  a  CW  laser  at  377  nm,  which  would  also  couple 
the  bound  state  with  the  dissociative  state.  The  Paladin  had  a  clean  Gaussian  profile, 
while  the  CW  laser  had  the  typical  elliptical  shape  of  diode  lasers  [183].  We  noticed 
that  at  some  point  the  UV  would  detach  electrons  from  the  surface  of  the  electrodes, 
which  would  in  turn  dissociate  molecule.  We  observed  in  the  spectra  very  fast  dis¬ 
sociation  rates  as  shown  in  Fig.  37.  The  dissociation  rate  with  the  UV  only  were 
faster  than  the  expected  dissociation  rate,  which  were  observed  with  the  paladin.  We 
attributed  to  electrons  stripped  from  calcium  that  coated  the  surface  of  the  electrodes 
in  the  trap.  The  work  function  of  calcium  is  2.9  eV,  while  the  energy  of  the  UV  laser 


is  3.2  eV,  which  means  the  laser  had  the  energy  required  to  strip  the  electrons.  The 
molecular  ion  could  then  be  dissociated  by  the  combination  with  the  electron  in  one  of 
the  two  ways.  Either  through  dissociative  excitation  (40CaH+  +  e“  — )-CaH*— >Ca++ 
H+  e“)  or  ion  pair  formation  (40CaH+  +  e“  — >-CaH*— >-Ca++  H~)  [214],  These  two 
processes  are  normally  endothermic,  however,  since  the  electrons  become  accelerated 
by  the  trapping  electric  field,  these  processes  can  be  observed.  We  did  not  conduct  a 
study  to  determine  the  actual  dissociative  mechanism  in  our  trap. 

To  get  around  this  problem,  we  opened  and  cleaned  the  trap  to  get  rid  of  the 
calcium  with  warm  sulfuric  acid.  Furthermore,  we  cleaned  up  the  mode  377  nm  laser 
to  minimize  scattering  of  the  laser  from  the  trap  electrodes.  This  reduced  the  number 
of  electron  bombardment  events  on  the  dissociation  rates.  An  electron  bombardment 
event  was  when  the  dissociation  time  was  less  ten  seconds,  which  was  faster  than  the 
expected  dissociation  time.  During  each  experiment,  we  checked  for  electron  bom¬ 
bardment  by  running  the  dissociation  experiments  with  the  UV  only  and  observing 
the  dissociation  times.  After  about  six  experiments  that  are  consistent  with  the  back¬ 
ground  dissociation  rate,  we  would  begin  REMPD  experiments  as  shown  in  Fig.  37. 
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Figure  37:  Figure  showing  the  effects  of  electron  recombination  on  the  dissociation 
rate.  Electron  recombination  resulted  in  dissociation  rates  that  were  significantly 
faster  than  the  observed  background  rates  and  significantly  higher  than  two-photon 
dissociation  rates  (Fig.  34).  The  processes  to  get  around  electron  recombination 
resulted  in  significant  lower  dissociation  rate  as  shown  in  the  figure. 


CHAPTER  VI 


FUTURE  DIRECTIONS  AND  GROUND  STATE 
COOLING  OF  A  MOLECULAR  ION 

6. 1  Introduction 

Resonance  enhanced  multiphoton  dissociation  is  a  great  tool  for  initial  survey  of 
molecular  ion  spectra  because  of  the  relative  ease  of  implementation.  However,  the 
limit  of  REMPD  is  the  linewidth  of  the  laser.  To  get  the  kind  of  precision  needed  for 
probing  fundamental  physics  (~1CT15),  higher  precision  measurements  are  needed. 
Such  a  technique  was  demonstrated  by  the  Wineland  group  at  the  National  Institute 
of  Standards  and  Technology  (NIST)  when  they  measured  the  clock  transition  of  an 
Al+  ion  co-trapped  with  Be+  [31].  The  method  called  quantum  logic  spectroscopy 
(QLS)  coupled  the  motion  of  the  ions  in  trap  with  the  internal  excitation  of  the 
ions  [31,  51,  212,  215,  216].  Thus,  the  clock  transition  of  Al+  was  measured  to  a  high 
accuracy  (1CT17)  without  the  need  to  directly  laser-cool  the  ion.  Such  a  technique  is 
congenial  for  spectroscopy  of  molecular  ions  because  most  molecular  ions  cannot  be 
directly  laser-cooled. 

However,  the  implementation  of  QLS  is  rather  difficult  as  it  requires  sophisticated 
experimental  control.  An  alternative  to  QLS  is  to  observe  the  heating  of  a  trapped 
crystal  caused  by  internal  excitation  and  spontaneous  emission  of  a  sympathetically- 
cooled  ion.  In  our  lab,  we  implemented  such  a  technique,  which  we  called  sympathetic 
heating  spectroscopy  (SHS)  [188].  Briefly,  two  ions  are  trapped;  one  is  directly  laser- 
cooled  (logic  ion)  while  the  other  is  sympathetically-cooled  (spectroscopy  ion).  A 
transition  on  the  spectroscopy  ion  is  probed  by  observing  the  heating  of  the  crystals 
as  a  function  of  the  detuning  from  that  transition.  While  the  probed  ion  scatters  a  few 


photons,  they  are  enough  to  cause  observable  heating  that  can  be  detected.  While 
SHS  works  well  for  electronic  transitions,  molecular  transition  scatter  even  fewer 
photons  because  of  the  availability  of  many  decay  channels.  Thus  a  more  sensitive 
technique  is  required.  For  this  purpose  we  proposed  observing  the  heating  of  sidebands 
of  ground  cooled  crystals  to  detect  molecular  transitions.  We  called  this  technique 
quantum  sympathetic  heating  spectroscopy  (QSHS).  Another  group  demonstrated 
this  technique  albeit  on  an  atomic  ion,  which  they  termed  photon  recoil  spectroscopy 
(PRS)  [217].  A  schematic  of  the  technique  is  shown  in  Fig.  38. 


Logic  Ion  Spec.  Ion 


(a)  Ground  motional  cooling 


Figure  38:  A  schematic  of  QSHS.  The  red  energy  spacing  lines  represent  the  motional 
state,  while  the  black  lines  represent  the  internal  state  of  the  molecule,  (a)  An  ion 
(logic  ion)-molecule  (spectroscopic  ion)  crystal  is  ground  state  laser-cooled,  (b)  A 
transition  in  the  spectroscopy  ion  is  probed  through  multiple  pulses,  which  causes 
emission  of  photons,  (c)  The  photons  cause  heating  of  the  sidebands  of  the  crystal, 
which  can  be  detected  (d)  as  a  function  of  the  detuning  of  the  wavelength. 

QSHS  works  by  cooling  the  crystal  to  the  ground  motional  state,  typically  the  axial 
mode.  A  transition  on  the  molecular  ion  is  probed  with  a  series  of  pulses.  Changes  in 
angular  momentum  of  the  spectroscopy  ion  is  transferred  to  the  logic  ion  through  the 
strong  Coulomb  attraction  in  the  ion-molecular  ion  pair.  Photon  scattering  at  the 
molecule  causes  heating  of  the  motional  modes  which  can  be  detected.  For  40CaH+, 
we  propose  to  frequency  double  the  TkSapph  laser  to  excite  the  electronic  transition 
of  40CaH+  of  the  21£+  •<—  11E+  electronic  transition.  The  u'—2  (393  nm),  v'— 3  (382 
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Figure  39:  The  proposed  experiment  to  measure  the  electronic  transition  21£+  t— 
11S+  of  40CaH+  taken  from  Ref.  [24], 

nm)  and  z/=4  (371  nm)  have  strong  Frank-Condon  overlap  with  the  ground  state. 
The  excited  molecule  should  fluoresce  around  616-686  nm.  This  photon  is  enough  to 
cause  significant  heating  in  the  sidebands,  which  can  be  monitored  as  a  function  of 
frequency  of  the  TkSapph  laser. 

A  prelude  to  QSHS  is  ground  state  cooling  of  the  ion-molecular  ion  pair,  which 
we  achieved  as  shown  in  the  next  section. 

6.2  Ground  state  cooling  of  a  molecular  ion 

Part  of  this  section  is  from,  Sympathetic  cooling  of  molecular  ion  motion  to  the  ground 
state 

R,  Rugango,  J.E  Goeders,  T.H.  Dixon,  J.M.  Gray,  N.  Khanyile,  G.  Shu,  R.  J.  Clark, 
K.  R.  Brown,  New  J.  Physics ,  17,  035009  (2015). 

6.2.1  Methods 

To  cool  the  motional  state  of  a  crystal,  narrow  transitions  on  a  metastable  state 
are  utilized  as  described  in  the  preceding  chapters.  Sideband  cooling  is  prefaced  by 
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Doppler  cooling  upon  which  the  ion  is  confined  in  a  space  much  smaller  than  the 
wavelength  of  the  transition.  A  laser  with  a  narrow  linewidth  can  be  tuned  to  the 
lower  motional  state.  After  each  excitation  to  the  excited  state,  the  vibrational  energy 
is  removed  as  shown  in  Fig.  40a.  Since  the  states  are  typically  long-lived,  to  ensure 
cooling  proceeds  faster  than  heating  processes,  the  long  lived  state  is  coupled  to  a 
fast  decaying  state  as  shown  in  Fig.  40b. 


Figure  40:  A  schematic  of  sideband  cooling,  (a)  An  excitation  to  an  upper  electronic 
state  removes  one  quanta  of  motional  energy  while  spontaneous  decay  results  in  no 
change  of  the  motional  state,  (b)  the  lifetime  of  the  D5/2  state  is  shortened  by  coupling 
it  to  the  P3/2  state,  which  spontaneously  decays  to  the  Si/2  state. 


For  calcium  ions,  the  S1/2  _D5/2  (729  nm)  transition  is  used  for  sideband  cooling 

experiments.  The  D5/2  state  has  a  relative  long  lifetime(l  sec),  the  cooling  process 
is  sped  up  by  coupling  to  the  P3/2  (854  nm)  state  for  fast  cooling.  Typically,  a 
magnetic  held  (~1  Gauss)  splits  the  degenerate  Zeeman  states  to  distinct  sublevels. 
The  transitions  for  sideband  cooling  are  then  chosen  as  the  <Si/2(m  =  —1/2) 

D5 /2 (rn  =  —5/2)  driven  at  729  nm  and  the  metastable  state  is  the  coupled  to  the 
P3/2(m  =  —3/2)  with  an  854  nm  laser.  The  P3/2(m  =  —3/2)  decays  primarily  back  to 
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the  S\/2{m  =  —1/2).  To  ensure  that  all  the  population  is  in  the  Si/2(m  =  —1/2)  state, 
spin  polarization  is  used  to  transfer  the  population  from  S\/2{m  =  +1/2)  state.  This  is 
achieved  by  driving  the  S\/2(rn  =  +1/2)  D5/2(m  =  —3/2)  then  again  depopulating 

it  with  an  854  nm  laser  to  the  to  the  P^/2{m  =  —1/2)  state  (Fig.  41).  This  process 
is  repeated  until  over  95%  of  the  population  is  transferred  to  the  S\/2{m  =  —1/2). 
Any  leakage  to  the  D3/2  state  is  depopulated  with  an  866  nm  that  is  always  kept  on 
during  the  entire  procedure. 
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(a)  Spin  polarization 
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(b)  Sideband  cooling  scheme 

Figure  41:  Spin  polarization  and  sideband  cooling,  (a)  shows  a  schematic  of  spin 
polarization  to  prepare  the  population  in  the  S\/2(m  =  —1/2)  state,  which  precedes 
sideband  cooling.  The  typical  sideband  cooling  scheme  is  shown  in  (b). 


Sideband  cooling  is  then  commenced  on  the  transition  stated  above.  The  electron 
shelving  technique  is  used  to  detect  the  state  of  the  ion.  Since  the  metastable  D  state 
is  non-fluorescent,  turning  on  the  Doppler  cooling  laser  can  reveal  whether  the  ion  is 
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the  D  state  or  the  S  state  by  the  observation  of  no  fluorescence  versus  fluorescence. 
A  typical  sequence  for  sideband  cooling  is  shown  in  Fig.  41b. 

6.2.2  Cooling  of  single  ion 


Prior  to  cooling  a  molecular  ion,  we  demonstrated  sideband  cooling  of  a  single  molec¬ 
ular  ion.  Briefly,  the  ion  was  trapped  in  a  trap  driven  at  14.426  MHz  resulting  in 
radial  secular  frequencies  of  1.419  MHz  (x-direction)  and  1.475  MHz  (y-direction). 
DC  voltage  confines  the  ion  in  the  axial  direction  with  a  secular  frequency  of  568 
KHz.  Upon  trapping,  the  ion  is  Doppler  cooled  and  sideband  cooled.  The  results  are 
shown  in  Fig.  42. 

To  determine  the  temperature  of  the  ion  after  sideband  cooling,  we  compare  the 
population  in  the  red  and  blue  sidebands.  Briefly,  after  sideband  cooling  the  average 
motional  state  (n)  is  in  a  Maxwell- Boltzmann  distribution  [129], 


(n)r 


Pn  \(n)  +  l)n+1 


(95) 


We  can  determine  the  population  in  the  red  sideband  by  applying  a  pulse  of  length 
t.  The  probability  that  the  ion  is  in  the  red  sideband  is, 

OQ  /  \Tl 

Prsb  =  ^2  /  /  \  j\n-|-l  sin  (^n,n-U/2)-  (96) 

n= 1  ''  ' 

For  the  blue  sideband,  the  population  is  given  by, 

oo  /  \n 

Pbsb  =  222  _|_  l)n+l  sin  (^n,n+U/2)-  (97) 

If  we  expand  the  sin2(f2f/2)  terms  around  0  and  then  use  the  Mclaurin  series  and 
compare  the  two  probabilities,  we  obtain  the  average  motional  quanta.  That  is, 


Prsb  (jl) 

Pbsb  (n)  +  1  ‘ 

Rearranging  the  above  equation  and  solving  for  (n),  we  find, 


(98) 
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Once  we  measure  the  population  distribution  of  sidebands  we  can  calculate  the  tem¬ 
perature  of  the  ion  by  the  simple  energy  relation, 


kBT  =  frw0({n)  +  ^), 
where  ujq  is  the  secular  frequency  of  the  trapped  ions. 


(100) 


Offset  from  carrier  (MHz) 


Figure  42:  Sideband  cooling  of  one  ion.  (a)  shows  population  in  the  red  sidebands  after 
sideband  cooling  while  (b)  shows  the  population  in  the  blue  sideband.  The  overall 
energy  of  the  system  is  found  by  comparing  the  two  peaks. 

After  Doppler  cooling,  the  sidebands  are  approximately  the  same  height.  The 
temperature  is  measured  to  be  T  =  0.75  mK  using  the  carrier  Rabi  oscillation  deco¬ 
herence  method  [218,  219],  which  is  close  to  the  Doppler  cooling  limit  of  T  =  0.53 
mK.  Sideband  cooling  is  performed  on  the  axial  mode  and  naxiai  =  0.1  is  typically 
achieved  as  measured  by  the  peak  height  comparison  method  [129] .  The  results  show 
that  the  motional  state  of  the  ion  after  sideband  cooling  is  found  to  be  n=0.1,  which 
corresponds  to  a  temperature  of  16.7±2.3  /rK. 

Sideband  cooling  of  a  molecular  ion  is  achieved  in  a  similar  fashion  as  sideband 
cooling  of  a  single  ion.  40CaH+  is  prepared  by  loading  two  ions  inside  the  trap  and 
Doppler  cooling  them.  IR  gas  is  introduced  into  the  chamber  via  a  leak  valve,  which 
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reacts  with  one  of  the  ions  in  the  trap.  This  is  shown  on  the  EMCCD  camera  by  a 
dark  spot  where  the  previously  fluorescing  ion  was  located.  The  non-reactant  ion  does 
not  change  position  indicative  of  the  presence  of  another  non-fluorescing  ion  in  the 
trap.  The  identity  of  the  molecular  ion  is  confirmed  using  the  method  described  in 
a  previous  experiment  [206].  The  sideband  cooling  scheme  proceeds  by  first  Doppler 
cooling  the  crystal  for  500  fx s  with  the  397  nm  laser  detuned  10  MHz  from  resonance 
and  then  continuously  exciting  with  the  729  nm  laser  alternating  between  the  red  first 
order  COM  and  BM  sidebands  for  6  ms  with  the  854  nm  laser  on.  Each  cooling  cycle 
is  preceded  by  a  spin  polarization  phase  on  the  S1/2  (m  =  +1/2) — >D5/2  (m  =  -3/2) 
transition  to  prepare  the  ion  in  the  S1/2  (m=  -1/2)  Zeeman  state.  Spin  polarization  is 
interleaved  with  the  cooling  and  repeated  every  100  /is.  After  cooling,  we  probe  the 
two-ion  crystal  with  the  same  technique  used  to  measure  motional  sidebands.  The 
average  motional  quanta  after  cooling  is  determined  through  the  ratio  of  the  heights 
of  the  red  and  blue  sidebands  [129]. 

The  ion  crystal  occasionally  suffers  from  collisions  with  residual  background  gas. 
Collisions  during  the  sideband  cooling  procedure  will  melt  the  crystal  and  result  in 
temperatures  far  above  the  Doppler  cooling  limit.  Cooling  lasers  close  to  resonance 
may  not  completely  bring  the  crystal  back  to  temperatures  near  the  Doppler  limit 
within  the  normal  cooling  period.  Collision  events  are  detected  by  observing  below 
normal  fluorescence  during  the  regular  Doppler  cooling  stage.  After  detection  of  a 
collision,  a  short  frequency  chirped  Doppler  cooling  pulse  is  applied  that  recrystallizes 
the  ion  chain  and  achieves  the  desired  initial  temperature  for  performing  sideband 
cooling.  We  discard  the  electron  shelving  data  recorded  coincident  with  the  collision 
event  from  out  final  data  set. 
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Figure  43:  Measurement  of  the  S1/2  (nr  =  -1/2) — )-D5/2  (nr  =  -5/2)  first-order  axial 
sidebands  for  the  40Ca+-  40CaH+  crystal  by  electron  shelving.  Red  and  blue  side¬ 
bands  are  fit  to  a  Gaussian  of  the  same  width  but  variable  amplitude.  Comparison 
of  sideband  heights  yields  an  average  mode  occupation  of  ucom  =  0.13  ±  0.03  for  the 
center  of  mass  mode  and  ubm  =  0.05  ±  0.02  for  the  breathing  mode. 

Our  molecular  ion  results  are  presented  in  Fig.  43.  After  sideband  cooling,  we 
observe  that  the  red  peak  height  is  greatly  suppressed  relative  to  the  blue  peak  height. 
Peak  height  comparison  reveals  nCOM  =  0.13  ±  0.03  and  ?Xbm  =  0.05  ±  0.02.  To 
determine  the  temperature,  we  match  the  expected  occupation  of  the  oscillator  as 
a  function  of  temperature  with  the  measured  occupation.  We  find  TCOm  =  12.47  ± 
0.03  /iK  and  Tbm  =  15.36  ±  0.01  /iK,  which  is  more  than  a  factor  of  3  below  the 
Doppler  cooling  limit.  The  presented  data  is  typical  and  similar  results  were  seen 
even  when  imperfect  compensation  resulted  in  a  shifted  breathing  mode  frequency 
[74,  206].  Ground  state  cooling  of  the  motion  in  three  dimensions  was  not  achieved 
due  to  an  unusually  high  radial  heating  rate.  We  were  only  able  to  cool  to  rTradiai  =  1 
corresponding  to  30%  population  in  the  ground  state.  This  could  be  due  to  residual 
micromotion  or  stronger  coupling  between  the  fluctuating  electric  fields  and  the  radial 
direction. 
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The  sideband  cooling  results  are  an  important  step  towards  the  implementation  of 
quantum  logic  spectroscopy  [31]  or  photon-recoil  spectroscopy  [217].  These  techniques 
are  limited  by  trap  heating  which  acts  as  unwanted  background  that  could  mask  the 
signal.  For  our  current  experiment,  we  measure  a  single  ion  background  heating  rate 
of  0.1  quanta/ms  comparable  with  other  experiments  in  similar  scale  traps  [220].  The 
heating  rate  of  the  atomic  and  molecular  ion  on  the  center  of  mass  mode  was  0.3 
quanta/ms.  If  this  mode  is  used  for  spectroscopy,  the  heating  rate  provides  an  idea 
of  how  quickly  the  molecular  ion  must  absorb  or  scatter  photons  to  have  a  detectable 
signal. 

6.3  Conclusions  and  future  outlook 

We  have  presented  the  first  measurements  on  the  vibrational  spectra  of  sympatheti¬ 
cally  laser-cooled  molecular  ion  40CaH+.  Guided  by  theoretical  calculations,  we  uti¬ 
lized  a  broadband  mode-locked  Ti:Sapph  laser  and  a  CW  laser  to  resonantly  dis¬ 
sociate  the  trapped  40CaH+.  Since  the  molecular  ion  had  an  internal  temperature 
equilibrated  with  black  body  radiation,  many  rotational  levels  are  thermally  popu¬ 
lated.  The  broadband  laser  effectively  addressed  all  these  rotational  levels  allowing 
only  the  vibrational  structure  to  be  probed.  Furthermore,  we  cooled  the  ion  down  to 
its  motional  ground  state  inside  the  trap,  which  makes  it  possible  to  perform  high 
precision  measurements.  We  propose  to  study  the  electronic  spectra  of  40CaH+  using 
the  heating  of  narrow  transitions  of  ground  state  cooled  molecular  ion,  a  method  we 
call  QSHS.  For  precision  measurements,  it  is  ideal  to  have  the  internal  population  of 
the  molecule  in  the  rovibrational  ground  state.  For  40CaH+,  the  vibrational  popula¬ 
tion  is  already  in  the  ground  state,  while  the  rotational  states  are  still  in  equilibrium 
with  black  body  radiation.  Cooling  the  rotational  motion  can  be  done  in  several 
ways.  Several  groups  have  utilized  optical  pumping  [66,  67],  which  has  a  disadvan¬ 
tage  of  requiring  lasers  that  are  deep  IR,  which  may  not  be  commercially  available. 


In  addition,  the  technique  requires  detailed  knowledge  of  the  transitions  involved.  An 
alternative  is  to  use  collisions  with  cold  gases  to  lower  the  internal  temperature.  This 
has  been  demonstrated  with  a  magneto-optical  trap  (MOT)  and  a  helium  cryogenic 
trap  [70,  73].  An  advantage  of  these  techniques  is  the  applicability  to  a  wide  range 
of  molecular  ions  without  the  need  for  detailed  spectroscopic  knowledge  of  those 
molecules.  In  the  Brown  lab,  the  construction  and  implementation  of  a  helium  buffer 
gas  cryogenic  cooling  has  been  demonstrated  on  a  chain  of  calcium  ions  [221],  The 
next  step  is  to  sympathetically  rotationally  cool  40CaH+  with  a  helium  buffer  gas. 
The  Brown  lab  is  also  interested  in  directly  laser  cooling  a  molecular  ion.  To  date, 
direct  laser  cooling  of  a  molecular  ion  has  not  been  demonstrated.  Few  candidates 
have  been  proposed  as  ideal  for  direct  laser  cooling  experiment,  BH+,  A1H+,  and 
SrO+  [44,  45].  For  this  experiment,  we  propose  to  cool  BH+  [13,  44,  222,  223].  Direct 
laser-cooling  of  molecules  requires  diagonal  Franck-Condon  factors,  which  would  re¬ 
quire  fewer  re-pump  laser  in  order  to  reach  the  translationally  cold  limit.  Direct  laser 
cooling  of  molecular  ions  reduces  unwanted  Stark  shifts  introduced  by  the  atomic 
cooling  lasers,  which  would  allow  for  higher  precision  measurements. 
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APPENDIX  A 


ANCILLARY  MATERIAL 

The  table  below  shows  the  some  physical  properties  of  the  transitions  in  40Ca+  ion. 
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Table  6:  Fundamental  physical  constants  used  in  this  thesis.  The  values  were  taken 
from  [430]. 


Fundamental  constant 

Symbol 

Value 

Speed  of  light  in  vacuum 

c 

299  792  458  m  s"1 

Planck  constant 

h 

6.62606957xl0-34  J  s 

Reduced  Planck  constant 

h 

1.054571726xl0-34  J  s 

Proton-electron  mass  ratio 

mp/me(n ) 

1836.15267245 

Fine-structure  constant 

a 

7.297  352  5698xRr3 

Electron  mass 

me 

9.10938291  xl0~34  kg 

Permittivity  of  free  space 

£o 

8.854187817 x  10-12  F/m 

Proton  mass 

mp 

1.672621777xl0-27  kg 

Elementary  charge 

e  or  q 

1.602176565X10”19  C 

Boltzmann  constant 

ks 

1.380  6488  xlO-23  J  K"1 
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